PORE    STUDIES    OF    CARBONATE     AGGREGATES 
BY   AN    EVAPORATION    TECHNIQUE 


SEPTEMBER    1968 
NO   22 


PURDUE    UNIVERSITY        \ 
LAFAYETTE  INDIANA  f 


Informational  Report 

PORE  STUDIES  OF  CARBONATE  AGGREGATES  BY  AN 
EVAPORATION  TECHNIQUE 


To:  J.  P.  McLaughlin,  Director 

Joint  Highway  Research  Project 

From:  H.  L.  Michael,  Associate  Director 
Joint  Highway  Research  Project 


September  12,  196B 
Project  No.: 
File  No.:  U-6 


The  attached  report  titled  "Tore  Studies  of  Carbonate  Aggregates 
by  an  Evaporation  Technique"  by  Harold  Reid  Brown  is  submitted  to  the 
Board  as  information.  Mr.  Brown  performed  the  research  under  the 
guidance  of  W.  L.  Dolch  of  our  staff  while  using  laboratory  facilities 
of  the  Portland  Cement  Association  and  while  employed  by  that 
organization. 

The  purpose  of  the  project  was  to  investigate  the  phenomena 
associated  with  the  evaporation  of  a  liquid  from  a  porous  solid,  to 
relate  these  to  pore  parameter  of  the  solid  and  to  apply  these 
techniques  to  carbonate  coarse  aggregates  in  an  effort  to  predict 
their  durability  in  concrete  under  freezing  and  thawing.  For  the 
latter  purpose,  laboratory  freezing  and  thawing  was  found  to  categorize 
the  aggregates  reasonably  well. 

The  report  is  presented  to  the  Board  as  information. 

Respectfully  submitted, 


Harold  L.  Michael  7S 

Associate  Director 


HLM:mz 

cc:  F.  L.  Ashbaucher 

W.  L.  Dolch 

W.  H.  Goetz 

W.  L.  Grecco 

G.  K.  Hallock 

M.  E.  Harr 


R.  H.  Harrell 
J.  A.  Havers 
V.  E.  Harvey 
G.  A.  Leonards 
F.  B.  Mendenhall 
R.  D.  Miles 
J.  C.  Oppenlander 


C.  F.  Scholer 
M.  B.  Scott 
W.  T.  Spencer 
H.  R.  J.  Walsh 
K.  B.  Woods 
E.  J.  Yoder 


Informational  Report 

PORE  STUDIES  OP  CARBONATE 
AGGREGATES  BY  AN  EVAPORATION  TECHNIQUE 


by 

Harold  Reid  Brown 
Graduate  Student 


Joint  Highway  Research  Project 
File:  k-6 


Purdue  University 
Lafayette,  Indiana 
September  12,  1968 


Digitized  by  the  Internet  Archive 

in  2011  with  funding  from 

LYRASIS  members  and  Sloan  Foundation;  Indiana  Department  of  Transportation 


http://www.archive.org/details/porestudiesofcarOObrow 


ii 


ACKNOWLEDGMENTS 

The  author  acknowledges  with  pleasure  and  gratitude  the  assistance 
and  guidance  given  by  his  major  professor,  Dr.  W.  L.  Dolch,  throughout 
this  study. 

Financial  assistance  for  this  work  was  made  possible  by  a  cooperative 
effort  between  the  Portland  Cement  Association  and  the  Joint  Highway 
Research  Project  of  the  School  of  Civil  Engineering,  Purdue  University, 
a  joint  activity  between  the  university  and  the  Indiana  State  Highway 
Department. 

The  author  is  grateful  to  Mr.  Paul  KLieger,  Manager,  Applied 
Research  Section,  Portland  Cement  Association,  for  his  timely  encourage- 
ment and  the  generous  support  of  his  staff  and  the  Association's  facilities, 

The  contributions  of  Dr.  Sidney  Diamond,  Dr.  C.  F.  Scholer,  and  Mr. 
Robert  Landgren  are  gratefully  acknowledged. 

Finally,  the  author  acknowledges,  with  deep  understanding  and 
appreciation  the  infinite  patience  and  encouragement  of  his  wife,  Evelyn, 
and  his  children. 


iii 


TABLE  OF  CONTENTS 

Page 

LIST  OF  TABLES v 

LIST  OF  FIGURES vi 

ABSTRACT ix 

INTRODUCTION  1 

Statement  of  Problem 1 

Objective  2 

Literature  Review  3 

Evaporation  Rates  9 

Approach  Used  in  This  Work 14 

EXPERIMENTAL  WORK 15 

Materials 15 

Sample  Preparation   17 

Apparatus  and  Procedure  19 

Absorption 19 

Evaporation  Rates   20 

Mercury  Porosimetry 22 

Tortuosity 25 

Tensile  Strength  27 

Freezing  and  Thawing 27 

Data 31 

General 31 

Evaporation  Rates   31 

Mercury  Porosimetry ♦ 39 

Freezing  and  Thawing 58 

DISCUSSION  .  . 63 

Theoretical  Considerations   63 

Constant -Rate  Period  63 

First  Falling  Rate  Period 63 

Second  Falling  Rate  Period  6k 

Permeability  .......  66 

Results .  67 

Discussion  of  Data  and  Results 72 

Evaporation  by  Air  Drying 72 

Evaporation  by  Vacuum  Drying  7^ 

Interpretation  and  Application  77 


iv 


SUMMARY  OF  RESULTS 
CONCLUSIONS  .  .  . 
BIBLIOGRAPHY  .  .  . 


Page 

90 
92 
93 


APPENDIX  A:      INFLUENCE  OF   ION  ADSORPTION  ON  TORTUOSITY 
MEASUREMENTS  FOR  MATERIALS  HAVING  HIGH 
INTERNAL  SURFACE  AREAS .        98 


APPENDIX  B:      THE  MATHEMATICAL  RELATIONSHIP  BETWEEN  THE  MASS- 
RATE   FUNCTION  AND  TIME-RATE   FUNCTION      


102 


APPENDIX  C 
APPENDIX  D 
APPENDIX  E 
APPENDIX  F 
APPENDIX  G 
APPENDIX  H 
APPENDIX  I 


NUCLEATION  RADII   103 

DRYING  DATA  FOR  VACUUM  EVAPORATION   105 

FREEZING  AND  THAWING  RESULTS   109 

GEOMETRY  CONSIDERATIONS  121 

FLOW  EQUATION ,  12U 

LIST  OF  SYMBOLS 127 

ASSUMPTIONS  OF  ANALYSIS  OF  AIR-DRYING  BY 

DIFFUSION 131 


VITA 135 


LIST  OF  TABLES 
Table  Page 

1.  Aggregate  Identification,  Source,  Geological  Origin 

and  Performance  ............  16 

2.  Porous  Materials,  Identification,  Description  and 

Properties , 1? 

3»   Coarse  Aggregate  Gradation 28 

k.     Fine  Aggregate  Gradation  .........  28 

5.  Concrete  Mix  Properties » »    29 

6.  Density,  Porosity,  Absorption,  Tensile  Strength,  Saturation 

and  Tortuosity 32 

?.  Insoluble  Residue  Contents  and  Qualitative  Identification 

by  X-ray  Diffraction 33 

8.  Slope,  "b",  for  Vacuum  Dried  Aggregates  and  Their 

Respective  Porosity,  Tortuosity  and  Area  .........    ^0 

9.  Calculated  Fore  Diameter  .........  68 

10.  Calculated  Permeabilities  69 

11.  Coefficients  of  Diffusion 71 

12.  Freezing  and  Thawing  Results 71 

13.  Saturation  and  Drying  Times  . 80 

1^.  Summary  of  Various  Effects  of  Freezing  Saturated 

Aggregates  ..........  83 

Appendix 
Table 

Dl.  Drying  Data  for  Vacuum  Evaporation 105 


vi 


LIST  OF  FIGURES 

Figure  Page 

1.  Typical  Drying  Rate  Curve 11 

2.  Air  Drying  Apparatus 21 

3.  Vacuum  Drying  Apparatus  ..................  23 

4.  Relation  between  Vacuum  Drying  Rate  and  Volume  Loss 

for  9-15,  I  A,  IS 34 

5.  Relation  between  Vacuum  Drying  Rate  and  Volume  Loss 

for  ML,   II 35 

6.  Relation  between  Vacuum  Drying  Rate  and  Volume  Loss 

for  47-2S,   IDU,  IA 36 

7.  Relation  between  Vacuum  Drying  Rate  and  Volume  Loss 

for  67-2S,  IDE,  MM 37 

8.  Relation  between  Vacuum  Drying  Rate  and  Volume  Loss 

for  MTS,  ML 38 

9.  Air  Drying  Rate  Curves  for  1-1S,  67-2S,  47-2S  and 

9-1S 41 

10.  Air  Drying  Rate  Curves  for  IS,  IA  and  II 42 

11.  Air  Drying  Rate  Curves  for  MTS,  ML  and  Brick 43 

12.  Air  Drying  Rate  Curves  for  IDE  and  IDU 44 

13.  Air  Drying  Rate  Curves  for  47-2S  Using  Toluene,  Isobutyl 
Alcohol  and  Sodium  Lauryl  Sulfate  Solution  45 

14.  Air  Drying  Rate  Curves  for  Fritted  Glass  Filters 46 

15.  Air  Drying  Rate  Curves  for  Cement  Fastes  and  Vycor  ....    47 

16.  Relation  between  Air  Drying  Rate  and  Volume  Loss  for 

1-1S,  67-2S,  IS 48 

17 •  Relation  between  Air  Drying  Rate  and  Volume  Loss  for 

IDE,  IDU,  II,  ML,  MTS 49 


vii 


Figure  Pag© 

18.  Relation  between  Air  Drying  Rate  and  Volume  Loss  for 

47-2S,  9-1S,  IA,  KM 50 

19.  Pore  Size  Distribution  for  1-1S,  4?- 23,  IDE 51 

20.  Pore  Size  Distribution  for  9-1S,  67-2S,  IDU 52 

21.  Pore  Size  Distribution  for  IS,  II,  MM 53 

22.  Pore  Size  Distribution  for  IA,  MTS,  ML 5^ 

23.  Fore  Size  Distribution  for  Brick 55 

24.  Pore  Size  Distribution  for  Oven  Dried  Cement  Pastes  ....  56 

25.  Fore  Size  Distribution  for  Vycor 5? 

26.  Freeze- Thaw  Expansion  of  1-1S,  4?-2S,  IA  59 

27.  Freeze- lhaw  Expansion  of  IDE,  IDU,  II 60 

28.  Freeze- Thaw  Expansion  of  ML,  MTS 61 

29.  Freeze- Thaw  Expansion  of  9-1S,  IS 62 

!.0.  Comparison  of  Years  of  Service  with  Expansion  Factor  ...  88 


•7 


Appendix 
Figure 

Al.  Absorbtivity  Curve  for  Vycor 100 

El.  Freezing  and  Thawing  Specimen,  MTS  ......  109 

E2.  Freezing  and  Thawing  Specimen,  ML 110 

E3.  Freezing  and  Thawing  Specimen,  9-1S  .....  Ill 

E4.  Freezing  and  Thawing  Specimen,  9-1S  (C)  .  .  .  . 112 

E5.  Freezing  and  Thawing  Specimen,  IS.. 113 

E6.  Relation  between  Relative  E  and  Freezing  and  Thawing 

of  MIS,  IDE,  IDU 114 

E7.  Relation  between  Relative  E  and  Freezing  and  Thawing 

of  1-LS,  9-1S,  II 115 

E8.  Relation  between  Relative  E  and  Freezing  and  Thawing 

of  47-2S,  IA,  IS 116 


vlii 


Appendix 
Figure  Page 

E9.  Relation  between  Weight  Change  and  Freezing  and  Thawing 

of  ML,  9-1S,  1*7-2S 117 

E10.  Relation  between  Weight  Change  and  Freezing  and  Thawing 

of  1-1S,  MTS 118 

Ell.   Relation  between  Weight  Change  and  Freezing  and  Thawing 

of  IDE,  II 119 

E12.   Relation  between  Weight  Change  and  Freezing  and  Thawing 

of  IS,  IDU,  IA 120 

Fl.  Schematic  of  Evaporation  from  a  Cylindrical  Specimen  ....  122 

II.  Typical  Examples  of  Air -Drying  Data  During  Second 

Falling  Rate  Period 131*- 


ABSTRACT 

Brown,  Reid  Harold.   Ph.D.,  Purdue  University,  June  1968.   Pore 
Studies  of  Carbonate  Aggregates  by  an  Evaporation  Technique.  Major 
Professor:  W.  L.  Dolch. 

The  purpose  of  this  study  was  to  investigate  the  phenomena  associated 
with  the  evaporation  of  a  liquid  from  a  porous  solid  and  to  relate  these 
to  pore  parameters  of  the  solid.   It  was  further  the  purpose  to  apply 
these  techniques  to  carbonate  rocks  used  as  concrete  coarse  aggregates 
in  an  effort  to  predict  their  durability  in  concrete  under  freezing  and 
thawing  conditions. 

A  variety  of  porous  media,  and  especially  a  suite  of  twelve  carbonate 
rocks  of  varying  field  histories  in  pavements,  were  used.   Small  samples 
were  saturated  under  vacuum  with  water.  The  water  was  removed  from  the 
pores  either  by  exposure  to  a  vacuum  or  to  a  dry  atmosphere  at  ambient 
pressure.  The  rates  of  drying  were  determined. 

The  porosities,  absorptions,  and  tortuosities  of  the  materials  were 
determined  by  conventional  techniques.  The  pore  size  distribution  curves 
were  determined  by  the  mercury  intrusion  method  using  an  apparatus  with 
a  maximum  pressure  of  5Q>000  psi. 

The  tensile  strengths  of  the  aggregates  were  determined  by  the 
cylinder  tensile  splitting  method.   The  aggregates  were  incorporated 
while  vacuum  saturated  in  air -entrained  concrete  prisms,  which  were  then 
subjected  to  cycles  of  freezing  in  air  and  thawing  in  water.  The  dynamic 


moduli  of  elasticity,  lengths,  and  weights  of  the  prisms  were  determined 
during  the  course  of  the  exposure. 

It  was  found  that  the  air-drying  process  gave  no  reliable  indication 
of  the  pore  parameters  of  the  solid,  and  specifically  the  rate  of  drying 
during  the  constant  rate  period  was  unrelated  to  pore  size.  An  empirical 
equation  for  this  process  was  developed  from  which  a  diffusion  coefficient 
was  calculated. 

For  the  vacuum-drying  process,  a  modified  slip-flow  equation  was 
developed  that  seems  to  describe  the  data  during  the  second  falling-rate 
period.   By  means  of  this  equation  average  pore  sizes  of  the  aggregates 
were  calculated  that  agree  reasonably  well  with  those  determined  by  direct 
measurement  of  permeability  and  with  the  logarithmic  mean  pore  sizes 
obtained  from  the  mercury  intrusion  results. 

The  laboratory  freezing  and  thawing  exposure  was  found  to  categorize 
the  aggregates  reasonably  well  in  terms  of  their  field  histories, 
especially  when  length  change  was  used  as  the  parameter  on  which  judgment 
was  based.  When  the  data  developed  here  were  combined  with  reasonably- 
assumed  parameters  the  behavior  of  all  but  one  of  the  aggregates  was 
found  to  be  explicable  in  terms  of  the  hypothesized  mechanisms  of  Verbeck 
and  Landgren. 


INTRODUCTION 

Statement  of  Problem 
The  study  of  the  performance  of  coarse  aggregate,  one  of  the 

predominant  constituents  of  portland  cement  concrete,  has  become  of 

increasing  importance.  Governmental  agencies  (l,  2,  3)  and  private 

interests  (  h   )  have  expressed  concern  about  the  depletion  of  known 

sources  of  satisfactory  aggregates,  especially  near  metropolitan  areas. 

The  magnitude  and  severity  of  the  problems  being  encountered  with  concrete 

aggregates  are  exemplified  by  the  priority  that  has  been  allotted  to 

aggregate  research  by  the  National  Cooperative  Highway  Research  Program.* 

Their  objective,  as  outlined,  in  one  of  their  recent  reports  (5  )  is: 

"Theories  have  been  established  to  explain  the 
phenomena  which  have  been  observed  in  studies  of  concrete 
subjected  to  alternate  cycles  of  freezing  and  thawing, 
and  test  methods  have  been  devised  for  experimental 
means  to  identify  aggregate  particles  which  undergo 
excessive  volume  change. 

"However,  none  of  these  studies  has  been  extended 
to  the  point  of  establishing  a  test  method,  or  methods, 
by  which  the  quality  of  aggregate  particles  can  be 
determined  on  a  positive,  reliable,  rapid,  and  economical 
basis.  This  need  has  always  been  urgent,  but  is 
particularly  so  today  because  of  the  rapid  depletion 
of  sources  of  supply  from  which  aggregates  of  'known' 
performance  have  been  obtained.  It  is  purely  conjectural 
as  to  how  truly  representative  of  aggregate  quality 


*  A  research  program  conducted  under  a  three-way  agreement 
entered  into  June  1962  by  the  National  Academy  of  Sciences, 
National  Research  Council,  the  American  Association  of  State 
Highway  Officials,  and  the  U.S.  Bureau  of  Public  Roads. 


this  'known'  performance  is  because  of  the  inadequacy  of 
currently  specified  experimental  measures  to  provide  data 
which  could  be  correlated  with  field  performance. 
Consequently,  such  adverse  effects  as  cracking,  spelling, 
and  pop-outs  have  undoubtedly  occurred  in  instances  where 
requirements  of  existing  soundness  tests  have  ostensibly 
been  met." 

Three  basic  approaches  have  been  used  to  identify  aggregate 
performance  in  portland  cement  concrete.  These  are:   (a)  evaluation  of 
field  or  in-service  performance,  (b)  laboratory  performance  tests,  and 
(c)  comprehensive  studies  of  the  material's  physical  properties.  Each 
approach  has  in  the  past  provided  valuable  insight  into  the  various 
aspects  of  confined  aggregate  performance  and  its  behavior  when  subjected 
to  the  forces  of  its  environment. 

This  research  has  been  primarily  concerned  with  the  pore  systems 
of  concrete  aggregates,  and  how  these  systems  are  related  to  D-line 
(deterioration  line)  cracking,  a  phenomenon  associated  with  freezing  and 
thawing  failure  of  the  concrete.  Circumstantial  evidence  points  to  a 
functional  relationship  between  an  aggregate's  performance  in  freezing- 
and- thawing  and  its  porosity,  pore  size  and  size  distribution,  tortuosity, 
and  pore  continuity.  Many  investigators  (  5»  6,  ?,  8,  9»  10,  11,  12,  13,  1*0 
have  applied  such  techniques  as  permeability,  absorption,  microscopy, 
mercury  porosiraetry,  and  vapor  adsorption  to  study  the  pore  structure  of 
porous  media. 

Objective 
In  an  effort  to  develop  a  test  method  that  fulfills  the  requirements 
of  simplicity,  rapidity,  economy,  and  reliability,  new  ways  were  sought 
to  measure  the  pore  characteristics  of  porous  systems.  In  1951  Messer  (15) 
reported  pore  characteristic  determinations  by  an  evaporation  method. 


Even  though  Messer's  work  had  several  unaccountable  results,  some  of  his 
basic  ideas  seemed  to  merit  further  consideration. 

In  earlier  studies  at  Purdue  University,  Dolch  (  16  )  measured 
evaporation  rates  on  a  few  aggregates,  but  no  attempt  was  made  to  analyze 
the  data.  A  preliminary  review  of  these  data  seem  to  show  some  functional 
relationship  between  the  evaporation  phenomena  and  the  structure  of  the 
pore  systems. 

The  works  of  Messer  and  Dolch  supported  with  some  preliminary  tests 
lead  to  the  following  postulation:  "There  is  some  functional  relation- 
ship between  the  evaporation  characteristic  of  porous  media  and  its  pore 
structure.*1  The  aim  of  this  study  was  to  identify  this  relationship, 
and  relate  the  findings  therefrom,  coupled  with  other  supporting  data, 
to  the  durability  field  performance  of  concrete  aggregate. 

Literature  Review 

This  review  of  literature  is  intended  to  cover  the  present  knowledge 
of  aggregate  as  a  porous  material  and  the  relationship  thereof  to  its 
freezing-and-thawing  durability  in  concrete. 

The  pore  structure  of  porous  media  and  the  properties  thereof  have 
been  of  prime  interest  to  the  geopetroleum  industry,  who  have  contributed 
significantly  to  many  of  the  fundamental  concepts  of  pore  theory,  hydro- 
statics, Darcy's  Law,  permeability  equations,  and  multiple  phase  flow. 
Muskat  (1?)  compiled  the  first  comprehensive  treatment  of  homogeneous  fluid 
flow  through  porous  media.  Maskat's  work  along  with  later  treatments  by 
Carman  (18  )  and  Scheidegger  (19 )  are  reoommended  as  reviews . 

Powers  et  al  (  20  )  investigated  the  permeability  of  portland  cement 
paste;  this  parameter  measures  the  resistance  to  flow  of  a  liquid  in 


porous  media.  Dolch  (21)  adapted  permeability  techniques  to  his  aggregate 
studies.  He  used  the  Kozeny-Carman  relationship,  which  is  based  on  an 
analogy  of  Poiseuille's  law  and  visoous  flow,  and  permeability,  to 
calculate  a  mean  pore  radius.   The  other  approach  used  by  Dolch  to 
determine  pore  sizes  by  permeability  was  referred  to  as  slip-flow.   Ihis 
concept  takes  into  account  the  slip  contribution  of  a  gas  in  its  flow 
through  a  capillary.  Absorptivity,  a  measure  of  the  rate  of  saturation 
by  capillarity,  and  surface  area  measurements  by  vapor  sorption  methods, 
were  also  used  by  Dolch  to  determine  pore  radii. 

Scheidegger  (19)  discussed  the  inadequacies  of  the  Kozeny-Carman 
equation.  One  of  his  primary  objections  was  the  tortuosity  factor, 
which  is  used  in  both  the  Kozeny-Carman  and  absorptivity  equations  to 
calculate  pore  sizes.  Dolch  (  21)  seemed  to  favor  the  absorptivity  method. 
Blanks  (  22)  used  the  capillary  suotion  phenomenon  on  rocks  and  on  concrete 
to  measure  pore  sizes.  With  the  materials  he  investigated,  a  diagnostic 
relationship  between  pore  characteristics  and  free zing-and- thawing 
durability  was  found. 

The  concept  of  mercury  porosiraetry  was  introduced  by  Washburn  (23) 
in  192?.  But,  it  was  not  until  19^5  that  Ritter  and  Drake  (2*+)  used 
this  approach  for  pore  studies.  Lemish  (12)  was  the  first  to  apply  this 
principle  to  concrete  aggregates  and  associated  materials.  Later,  he 
reported  with  Hiltrop  (13)  that  they  had  found  eight  typical  pore-size 
distribution  curves  among  the  carbonate  aggregates  they  had  studied. 
Their  conclusion  was  that  "Pore-size  distribution  alone  does  not  appear 
to  be  related  olosely  with  either  freezing-and-thawing  or  service  record" 
of  the  rocks  in  their  study,  and  that,  however,  some  relationship  does 
exist,  even  though  it  was  not  obvious,  and  further  work  was  in  order. 


The  absoissa  of  their  distribution  curve  was  plotted  with  a  linear  scale 
covering  a  range  from  zero  to  ten  microns.  With  most  of  the  aggregates 
they  studied,  the  majority  of  pores  were  smaller  than  one  micron.  What 
this  means  is  that  the  resolution  of  unique  and  distinguishable  features 
of  small  pores  on  the  distribution  curve  is  lost;  these  are  the  pores 
whose  influence  in  aggregate  performance  is  most  likely  to  be  significant. 

Two  other  studies  (  14,  25  )  have  recently  been  completed  using 
this  principle.  Sherwood  (  14  )  examined  fifteen  carbonate  aggregates. 
By  plotting  the  distribution  function  on  a  linear  scale  and  the  pore 
diameter  on  a  logarithmic  scale  he  was  able  to  distinguish  the  features 
of  the  pore  structure  for  small  pores  as  well  as  for  the  larger  ones. 
These  fifteen  carbonate  aggregates  were  grouped  into  four  classifications, 
according  to  their  pore  sizes.  Meininger  (25)  also  reported  similar 
pore  size  distribution  curves  on  carbonate  aggregates.  Neither  of  these 
reports  showed  significant  conclusions  between  pore-size  distribution 
and  durability  in  freezing  and  thawing. 

As  early  as  1936  (26)  petrographic  examinations  were  used  as  an  aid 
to  investigate,  select  and  control  concrete  aggregates.   This  method  of 
aggregate  analysis  is  widely  used  today  (  10,  26   )  and  contributes  much 
to  the  apperception  of  the  material's  properties,  but  is  limited  to  what 
is  observable.  With  a  light  microscope  it  is  difficult  clearly  to  resolve 
an  image  smaller  than  about  one  micron.  Sweet  (2?)  concluded  that 
aggregates  with  a  comparatively  large  void  volume  of  pores  smaller  than 
5  microns  in  diameter  were  likely  to  be  non-durable  in  freezing  exposure. 
If  the  voids  ratio  of  small  pores  was  10  percent  or  greater,  unsatisfac- 
tory performance  was  demonstrated,  while  those  materials  whose  voids 
ratio  of  small  pores  was  smaller  than  6  percent  showed  good  performance. 


However,  later  investigations  (  13»  ^»  °  )  showed  no  basis  for  this 
critical  index.  The  matter  is  discussed  further  by  Dolch  ( 6  ). 

Porosity,  e  ,  is  defined  as  the  ratio  of  the  pore  volume  to  the 
bulk  volume.  This  property  of  a  porous  material  has  been  shown  (8,  27) 
useful  in  evaluating  aggregate  durability.  Lewis  and  Dolch  ( 7  )  dis- 
cussed porosity  and  its  relationship  to  the  concrete  system.  Dolch  ( 6  ), 
in  a  subsequent  report,  classified  the  properties  influenced  by  porosity 
into  two  general  categories.  First  the  presence  of  pores  decreases  the 
total  solid  volume,  which  influences  properties  such  as  strength.  The 
second  property,  that  of  allowing  ingress  and  retention  of  water,  is 
perhaps  more  important  to  the  problem  of  durability  under  freezing. 

It  can  thus  be  seen  that  several  methods  are  available  that  will 
aid  in  describing  the  pore  system  of  porous  materials,  and  that  many 
attempts  have  been  made  to  relate  these  findings  to  the  performance  of 
aggregate.  However,  at  the  present  stage  of  development,  no  one  method, 
nor  even  a  combination  of  methods,  is  able  clearly  to  define  the  future 
durability  performance  of  an  aggregate. 

What  are  the  reasons  behind  this  apparent  failure  to  be  able  to 
answer  the  question:  How  should  this  aggregate  perform?  All  investi- 
gators have  recognized  the  fact  that  aggregates  are  inherently  composed 
of  a  complex  mixture  of  inhomogeneous  materials  and  this,  in  itself, 
complicates  the  investigation  process.  One  is  never  quite  sure  exactly 
what  he  may  have,  which  uncertainty  means  extensive  retesting  and 
statistical  treatments.  Whether  one,  two  or  even  several  parameters 
will  characterize  an  aggregate  is  still  hypothetical.  However,  what  are 
the  other  alternatives?  In  general,  two  other  approaches  have  been  used. 


A  significant  amount  of  work  (3,  27,  28,  29,  30,  31,  32,  33) has 
been  performed  with  concrete  in  a  laboratory  accelerated  freezing-and- 
thawing  environment.  This  test  method  is  widely  used  as  an  acceptance- 
rejection  test.  In  general,  if  a  specimen  shows  adequate  durability  in 
such  a  test,  its  performance  in  the  field  is  assumed  to  be  satisfactory. 
Of  course,  there  are  exceptions,  and  several  conflicts  remain;  these  are 
the  answers  we  seek.  o 

Some  investigators  (  29,  3^  )  have  contended  that  an  accelerated 
freezing-and- thawing  test  is  too  severe,  and  for  this  reason  many  mate- 
rials are  rejected  that  would  have  performed  well  in  the  field. 
Powers  (29)  proposed  an  empirical  test  that  may  be  more  relevant  to  field 
conditions  than  accelerated  tests.  Tremper  and  Spellman's  (  3^  ) 
results,  using  the  method  of  Powers  were  employed  as  a  basis  for  allowing 
questionable  aggregates  to  be  incorporated  into  a  California  highway 
(  y\   ),  and  these  aggregates  have  thus  far  proven  to  be  satisfactory. 
Wills  (35)  concluded  that  permanent  length  change  was  a  sensitive 
indicator  of  frost  damage  with  dilation  measurements  and  relative 
modulus  of  elasticity,  JS,  being  less  sensitive. 

It  is  the  consensus  of  the  majority  of  investigators  (2,  4,  29, 
30,  36  )  that  destructive  stress  emanates  from  the  cyclic  freezing- 
and-  thawing  of  water  within  the  aggregate-cement  pore  system  of  concrete, 
the  two  currently-accepted  mechanisms  of  this  action  are,  briefly,  as 
follows.  When  water  is  frozen  in  these  void  systems,  dilation  may  occur 
by  either  hydraulic  pressure,  growth  of  ice  lenses,  or  by  a  combination 
of  both  processes.  Hydraulic  pressures  are  generated  by  movement  of 
unfrozen  water  during  the  freezing  process  if  the  void  system  is  satu- 
rated above  the  critical  value  of  about  90  percent.  The  exact  magnitude 
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of  the  pressure  generated  depends  on  the  permeability,  air  content,  and 
degree  of  saturation  of  the  paste,  and  the  permeability,  porosity,  degree 
of  saturation,  and  size  of  the  aggregates.  Ice  lensing  is  a  process  in 
which  unfrozen  water  will  move  to  ice  within  the  porous  system  because 
of  a  free  energy  difference  between  the  two  states.  This  phenomenon  is 
common  in  soils  and  is  capable  of  developing  significant  pressures. 

Dunn  and  Hudec  (37)  have  recently  hypothesized  another  possible 
origin  for  destructive  stresses  in  a  concrete  system  when  subjected  to 
alternate  cycles  of  cooling  and  warming.  They  observed  by  cold  ETA 
methods  that  less  water  was  frozen  in  those  aggregates  classified  as 
poor  than  in  those  that  were  classified  as  good.  By  relating  these 
findings  to  the  adsorption  characteristics  of  the  aggregate  they 
concluded  that  the  destructive  stresses  were  originating  from  cycles  of 
wetting  and  drying  and  from  thermal  fluctuations.  The  aggregates  they 
were  investigating  had  high  argillaceous  contents.  This  approach  has  not 
so  far  met  with  wide  acceptance  (38) . 

Larson  et  al  (l)  concluded  that  the  most  promising  methods  to  pre- 
dict the  freezing-and- thawing  resistance  of  aggregates  were  by  porosity, 
permeability,  capillarity,  and  mercury  porosimetry.  Of  these  methods 
Larson  et  al  favored  mercury  intrusion,  from  which  pore  entrance 
diameters  and  their  distribution  can  be  calculated,  for  investigating 
the  pore  systems  of  concrete  aggregates. 

Verbeek  and  Landgren  (  k   )  stressed  the  importance  of  moisture 
migration  to  and  from  the  concrete  system.  The 'physical  characteristic 
of  the  aggregate,  its  porosity  and  pore- size  distribution,  and  the 
permeability  and  thickness  of  mortar  covering  the  aggregate  are  important 
factors  that  influence  time  requirements  for  saturation.  The  time  factor 


varies  with  seasonal  influences  -  drying  occurs  during  the  summer 
months  and  absorption  of  water  during  the  winter  or  wet  seasons.  An 
aggregate's  performance  may  be  entirely  satisfactory  if  the  environmental 
factors  are  favorable  to  drying  and  if  the  pore-systems  of  aggregate, 
paste,  or  concrete  are  respondent  to  drying;  whereas  if  the  environment 
had  been  unfavorable  to  drying  the  aggregate  or  concrete  may  have  failed. 

This  concept  offered  by  Verbeck  and  Landgren  suggested  a  possible 
new  approach  to  identifying  the  potential  performance  of  concrete 
aggregates.  Previous  studies  have  primarily  involved  themselves  with  how 
the  pore-system  acquires  its  moisture,  and  how  this  moisture  is  retained 
and  transported  through  the  system.  But  now  it  seems  that  an  under- 
standing of  how  the  moisture  leaves  a  system  merits  consideration.  One 
of  the  predominant  mechanisms,  by  which  moisture  is  removed  from  concrete, 
is  evaporation.  Undoubtedly,  the  process  of  removal  from  confined 
aggregates  and  from  unconfined  aggregates  is  different,  but  a  qualitative 
assessment  of  drying  from  unconfined  systems  will  certainly  aid  in  the 
overall  evaluation  of  moisture  transfer. 

Evaporation  Rates 

The  science  of  drying  has  been  extensively  investigated  and  reported, 
primarily  in  the  chemical  engineering  literature.   Industries  such  as 
food  processing,  mining,  and  textiles  have  contributed  much  to  drying 
technology.  The  soil  physicist  has  an  interest  in  the  drying  process. 
The  latest  theoretical  views  of  the  drying  process  are  summarized  by 
Tsao  and  Wheelock  (39)- 

Two  concepts  have  been  advanced  as  the  mechanism  by  which  fluids 
are  transferred  to  the  surface  for  evaporation.  One  is  that  the  liquid 


10 

flows  from  the  interior  to  the  outside  evaporating  surface  by  capillarity 
and  is  activated  by  a  pressure  difference  controlled  by  the  size  of  the 
pores.  Thus  the  rate  of  evaporation  should  be  a  function  of  the  pore  size 
(frictional  forces)  and  the  relative  humidity  (vapor  pressure  difference) 
between  the  interior  and  the  exterior  of  the  porous  body.  The  second 
mechanism  (40)  used  to  explain  moisture  transfer  is  diffusion.  The 
process  assumes  that  the  fluid's  surface  retreats  inwardly  as  drying 
proceeds,  and  that  the  moisture  moves  outwardly  as  a  vapor.  Fick's  law 
is  the  basic  equation  for  this  process. 

Derjaguin  (41)  showed  that  evaporation  rates  from  large  single 
capillaries  were  lower  than  those  from  smaller  ones,  which  is  contrary 
to  the  Kelvin  formula.  He  showed,  theoretically  and  experimentally,  that 
a  combination  of  film  flow  along  the  tube  surface  and  diffusion  through 
the  bulk  pore  space  explained  this  apparent  anomaly.   However,  Derjaguin 
did  not  consider  the  theory  applicable  to  a  multi-size  system  similar  to 
the  pore  system  of  concrete  or  aggregate  material.  He  made  reference  to 
Rallaire's  (42)  work  for  multi-size  systems. 

Drying  data  are  frequently  plotted  as  the  rate  of  drying  of  the  porous 
body  against  moisture  content,  Fig.  1.  The  plotted  curve  will  generally 
consist  of  three  regions.  The  constant  rate  period,  AB,  is  an  interval 
in  which  the  surface  of  the  porous  medium  remains  sufficiently  wet  to 
simulate  evaporation  from  an  equivalent  body  of  water  and  is  essentially 
Independent  of  the  porous  body's  properties.  This  constant  rate  period 
is  followed  by  a  period  of  evaporation  from  a  saturated  surface  of  gradu- 
ally decreasing  area,  and  diffusion  of  the  resulting  vapor  through  the 
pores  to  the  outside.  This  region  is  referred  to  as  the  first  falling-rate 
period,  BC.  The  shape  of  the  second  falling  rate  period,  CD,  is  governed 
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by  internal  diffusion  and  depends  upon  the  geometric  shape  of  the  specimen. 
The  point  of  inflection  between  the  constant  rate  period  and  the  first 
falling  rate  period  is  called  the  critical  point. 

The  process  of  moisture  movement  to  the  drying  surface  has  been  the 
subject  of  much  debate.  Ceagiske  and  Hougen  (40)  established  that  mois- 
ture is  moved  by  a  capillary  process,  which  finding  was  in  direct  oppo- 
sition to  earlier  theories  (39)  of  liquid  diffusion.  They  employed  a 
suction  potential  method  to  formulate  their  hypothesis.  This  concept  was 
used  later  by  Blanks  (22)  on  concrete  and  aggregates.  This  mechanism  is 
of  special  interest  to  the  soil  physicist  and  has  been  reviewed  by  Keen 
(*3). 
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Hallalre  (42.),  using  a  porous  soil  model,  identified  two  transfer 
potentials,  or  two  types  of  diffusion  (micro  and  macro).  He  developed  a 
theoretical  model  for  an  evaporation  mechanism  for  a  porous  bead  pack,  a 
model  similar  to  a  soil  structure.  His  work  attempted  to  explain  the 
anomalies  between  capillarity  and  diffusion  concepts.  His  model  had  two 
Independent  systems,  micro  and  macro,  whereby  he  was  able  to  explain  the 
apparent  abnormality  of  moisture  movement  from  low  moisture  gradients  to 
higher  gradients. 

Messer  (15);  in  an  effort  to  determine  interstitial  water  contents 
for  petroleum  bearing  rocks,  used  the  principle  of  evaporation.   He  dis- 
cussed in  some  detail  the  relationship  between  evaporation  rates  and 
pore  sizes,  and  used  a  theoretical  treatment  based  on  the  Kelvin  equation, 

P    eXp      r^RT  V' 

o 

where:   P/P  is  the  relative  vapor  pressure  over  a  liquid  having  a  sur- 
face tension  X,  molecular  weight  M,  density  p   ,  temperature  T,  and  a 
contact  angle  of  8  in  a  pore  of  radius  r,  and  R  is  the  gas  constant. 

Bartlett  (44)  suggested  a  relationship  between  evaporation  character- 
istics and  pore  size.  He  made  the  following  conclusions:  The  major 
critical  point  occurs  when  the  pores  of  modal  and  greater  sizes  have  been 
emptied  by  capillarity.  This  point,  and  the  effective  pore  size  distri- 
bution are  dependent  on  the  nominal  rate  of  drying.  These  conclusions, 
however,  were  based  partially  on  capillary  suction  tests  and  deductive 
logic  and  not  on  measurements  of  direct  evaporation. 

The  evaporation  process  is  complex  in  nature,  especially  in  a  porous 
system  where  several  parameters  are  involved.  Since  In  this  work  the 
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primary  interest  Is  to  identify  the  pore  properties  that  influence 
evaporation,  the  constant  rate  period  and  the  first  falling  rate  were 
not  used  in  the  analysis  of  the  data.  The  analysis  was  based  only  on  the 
second  falling-rate  period.  The  reasons  for  this  selection  are  given  in 
detail  in  the  Discussion  section. 
B'     The  differential  equation  of  diffusion  is 

^-DV^T  (2) 

0  t 

where  t  is  time,  C  the  concentration  of  evaporable  moisture,  D  the 
diffusion  coefficient  andV  the  Laplacian  operator.  Several  solutions 
to  equation  2  are  available  (  ^5,46).  By  applying  the  boundary  conditions 
of  a  uniform  distribution  of  moisture  having  a  concentration  falling  to 
zero  at  the  surface,  the  average  moisture  distribution  W,  for  a  long 
cylinder  of  length  L,  opened  at  one  end,  can  be  expressed  by  the  following 

series : 

n 
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Bancox  (4?)  by  ignoring  higher  harmonics  differentiated  equation  3 
with  respect  to  t. 


* ~   -  GDW  (5) 


or 

dt 


where  G  -  (|-)2 


Thus  if  the  diffusion  process  governs  and  harmonics  are  neglected,  the 

second  falling  rate  period  can  be  expressed  by  equation  5-  A  linear 

dW 
relationship  exists  between  the  rate,  —  and  W  if  D  is  independent  of  W. 


Ik 

This  is  then  the  means  whereby  the  diffusion  coefficient  can  be  determined 
from  the  data  of  the  seoond  falling-rate  period. 

Approach  Used  in  This  Work 

A  variety  of  porous  materials  were  selected  for  this  work  in  addition 
to  concrete  aggregates.  Vycor*  thirsty  glass,  fritted  glass,  sintered 
brass  filters,  and  three  cement  pastes  of  different  water-cement  ratios 
were  included  in  the  investigation.  Twelve  carbonate  aggregates, 
selected  primarily  on  the  basis  of  their  field  performances  and  test 
histories,  were  obtained  from  several  midwestern  states.  This  aggregate 
suite  provided  a  broad  spectrum  of  materials  the  past  performance  of 
which  in  concrete  pavements  ranged  from  good  to  bad.  The  materials  that 
were  classified  as  bad  had  shown  severe  lUcracking  in  concrete  pavements. 

True  density,  porosity,  tortuosity,  tensile  strength,  and  absorption 
were  determined.  Mean  pore  size  and  pore  size-distribution  were  measured 
using  the  mercury  intrusion  principle.  The  results  from  these  data  were 
compared  with  the  results  determined  by  evaporation. 

An  accelerated  freeze  in  air  and  thaw  in  water  test  for  concrete 
prisms  was  used  to  investigate  the  influence  of  these  coarse  aggregates 
on  the  durability  of  concrete,  and  to  assign  a  quantitative  durability 
rating  to  the  aggregate.  These  findings  were  related  to  the  physical 
properties  of  the  aggregates. 


*  Corning  Glass  Works,  Corning,  N.Y. 
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EXPERIMENTAL  WORK 


Materials 


In  the  past,  extensive  investigations  have  been  made  on  several 
Indiana  aggregates.  Four  of  these  were  selected  for  this  study.  Their 
geological  origin,  chemical  and  mineralogical  composition,  porosity 
absorptivity,  tortuosity,  permeability,  specific  surface,  and  field 
performance  were  reported  by  Dolch  (8).  The  quarry  data  were  reported  by 
Sweet  (48).  Portions  of  these  specimens  that  were  sampled  for  earlier 
studies  were  available  and  were  incorporated  in  this  study.  By  using 
these  samples,  the  problem  of  selecting  materials  representative  of  those 
from  former  studies  was  minimized,  and  thereby  the  results  of  both  studies 
can  be  compared  with  more  confidence. 

A  suite  of  additional  aggregates  was  selected  from  Iowa.  The 
physical,  chemical  and  engineering  properties  were  reported  by  Welp  and 
DeXoung  (49),  With  the  exception  of  the  Decorah  stone,  these  samples 
were  taken  from  ledges  currently  being  used  for  concrete  pavements.  The 
Decorah  stone  that  is  presently  used  shows  signs  of  deterioration  in 
pavement  after  about  20  years,  rather  than  the  age  earlier  reported  (49) 
of  10  to  15  years.  Three  aggregates  were  selected  from  other  sources. 
The  general  description  of  all  aggregates  incorporated  into  this  study  is 
summarized  in  Table  1. 

To  enlarge  the  spectrum  of  pore  size  and  structure  systems  of 
porous  media,  (carbonate  rooks  have  a  range  of  pore  sizes  usually  between 
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Table  1 

Aggregate  Identification,  Source,  Geological 
Origin  and  Performance 


Identt  f i- 

Geological 
Origin 

Performance 

cat ion- 
Cede 

City  and  State 

History 
in  Pavements 

1-1S 

Linn  Grove, 
Indiana 

Silurian 
(Doloraitic 
Limestone) 

Good 

67-2S 

Greenoastle , 
Indiana 

Mississippian 
(Limestone) 

Good 

k7-2S 

Bedford, 
Indiana 

Mississippian 
(Limestone) 

Poor 

9-1S 

Logansport, 
Indiana 

Silurian 
(Doloraitic 
Limestone) 

Poor 

MM 

Mankato, 
Minnesota 

(Limestone) 

Poor 
(10  years) 

ML 

Kansas  City, 
Kansas 

Argentine 
(Limestone) 

Poor 
(3  to  4  years) 

MTS 

Kansas  City 
Kansas 

Mississippian 
St.  Genevieve 

Poor 
(8  to  10  years 

(Limestone) 


IS 
II 

IDU 

IDE 
IA 


Menlo, 
Iowa 

Iowa  City, 


Dubuque, 
Iowa 


Decorah, 
Iowa 

Aid  en 


Argentine 
(Limestone) 


Poor 
(2  years) 


Devonian  Good 

Cedar  Valley  (40  +  years) 
(Limestone) 

Ordovician  Good 

Galena  (Up  to  35  years) 
(Dolomite) 


Ordovician 
Galena 


Poor 
(14  years) 


Mississippian      Good 
St.  Genevieve   (40  +  years) 
(Limestone) 
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0.01  and  10  microns)  other  materials  were  selected  for  study.  Vycor*,  a 
manufactured  porous  glass,  has  several  unusual  properties,  which  provide 
interesting  information  for  a  study  of  this  nature.  Other  porous 
materials  used  in  this  study  were  pyrex  fritted  glass  of  several  porosity 
grades,  sintered  brass,  hardened  portland  cement  paste,  and  ordinary 
brick.  Some  properties  of  these  materials  and  a  code  for  their  reference 
are  listed  in  Table  2. 

Sample  Preparation 

The  stones,  when  removed  from  the  quarries  were  large  and  irregular 
in  shape.  These  pieces  were  cut  parallel  to  the  bedding  plane  (the 
plane  was  identified  by  the  sampler  or  in  a  few  instances  the  bedding 
plane  was  assumed  by  visual  examination)  with  an  18- in.  Felker**  diamond 
saw  into  approximately  3- in.  slabs.  This  size  provides  for  easy  mounting 
of  the  specimen  in  a  modified  portable  Clipper***  drill  press  for  coring. 
The  Felker  core  bit  had  a  0.81 3- in.  inside  diameter. 

This  coring  bit,  which  is  used  primarily  for  cutting  glass  lenses, 
was  preferred  to  those  designed  for  heavy  duty  concrete  or  rock  coring. 
The  surface  of  the  cylindrical  wall  of  the  core  is  out  with  an  extremely 
smooth  surface  and  close  tolerances.  After  coring,  the  ends  of  the 
specimens  were  cut  off  perpendicular  to  the  cylindrical  wall  and  lapped 
on  a  glass  plate  to  a  final  finish  with  M-180  abrasive.  Each  specimen 
was  then  washed  for  five  minutes  in  an  ultra-sonic  washer****,  and  then 
oven-dried  at  105  C  for  2k  hr  before  being  stored  in  a  sealed  container. 


•*  Corning  Glass  Co.,  Corning,  N.Y. 
**  Felker  Mfg.  Co.,  Torrance,  Calif. 
***  Clipper  Mfg.  Co.,  Kansas  City,  Mo. 
****  Harris  Transducer  Corp.,  Woodbury,  Conn, 
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Table  2 

Porous  Materials,  Identification, 
Description  and  Properties 

Identification            Description  and  Properties 
Code  '  

Vycor  Porous  Glass 

Porosity  =  28$  2 

Internal  Surface  Area  =  250  cc  /gram 
Avg  Pore  Diameter     =  40  A 

Fire  clay  face  brick 

Fritted  glass  filter,  coarse 

Fritted  glass  filter,  medium 

Fritted  glass  filter,  fine 

Fritted  brass  filter,  coarse 

Portland  oement  paste,  hydrated  for 
11  years,  w/c  =0.35 

Portland  cement  paste,  hydrated  for 
11  years,  w/c  =0.65 

Portland  cement  paste,  hydrated  for 
11  years,  w/c  =  0.80 


Brick 

FGC 

FGM 

FGF 

FBC 

W/C  = 

0.35 

w/c  = 

0.65 

W/C  st 

0.80 
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The  material  used  for  mercury  intrusion  testing  was  selected  at 
random  from  the  prepared  samples.  They  were  crushed  to  a  size  that  passed 
a  No.  8  U.S.  Standard  sieve  and  was  retained  on  a  No.  16.  This  material 
was  rewashed  in  the  ultra-sonic  washer  for  another  five  minutes,  redried, 
and  bottled. 

All  samples  of  regular  geometrical  shape  were  measured  with  a 
micrometer  caliper  to  the  nearest  0.01  cm. 

Apparatus  and  Procedure 
Absorption 
Absorption  is  defined  as  the  ratio  of  the  weight  of  water  contained 
in  the  pores  of  a  material  fander  specified  conditions)  to  the  oven-dry 
weight  of  the  solid  material.  All  samples  were  evacuated  by  a  mechanical 
vacuum  pump  for  six  hours,  and  then,  while  under  vacuum,  were  completely 
submerged  with  the  fluid  under  study.  Then  the  vacuum  was  released  and 
the  sample  remained  submerged  at  atmospheric  pressure  for  Zh   hours.  For 
most  tests,  water,  which  was  distilled  and  de-aired  by  boiling  just  prior 
to  use,  was  the  saturant  fluid.  Others  included  toluene,  5  milliraolar 
solution  of  reagent  grade  and  specially  purified  sodium  lauryl  sulfate, 
and  isobutyl  alcohol.  This  technique  is  usually  referred  to  as  vacuum 
saturated,  and  is  considered  ( 6  )  to  saturate  a  porous  material  as 
completely  as  any  available  method.  The  sample,  after  removal  from  the 
fluid,  was  blotted  to  a  surface-dry  state,  a  point  where  the  surface 
appears  dull,  and  weighed  rapidly  on  an  analytical  balance.  The  retained 
fluid  is  referred  to  as  the  saturated  surface  dry  (SSD)  content. 
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Evaporation  Rates 

Two  drying  techniques  were  developed  to  investigate  the  evaporation 
characteristics  of  porous  media,  and  they  will  be  referred  to  as  the 
"Air  Dry"  and  "Vacuum  Dry"  methods. 

For  air  drying,  such  parameters  as  air  velocity,  temperature,  humidity, 
size,  shape,  and  time  had  to  be  controlled  and  varied  as  needed.   It  soon 
became  apparent  that  measurements  could  not  be  made  manually  because  of 
the  excessive  time  involved.  This  restriction  led  to  the  development  of 
an  automated  system. 

An  existing  thermo-balance  was  modified  to  measure  time-weight 
relationships  automatically.   This  equipment,  shown  by  a  photograph  in 
Figure  2  proved  to  be  completely  satisfactory  for  this  work.   The  sample 
was  suspended  in  the  stream  of  dry  air  coming  from  the  drying  train  and 
issuing  from  the  stack  above  the  balance.  The  weight  was  continuously 
recorded  on  the  recorder  and  was  accurate  to  0.1  mg.  The  final  result 
was  a  permanently  recorded  time-weight  curve  that  was  manually  differen- 
tiated to  obtain  evaporation  rates. 

The  sample  was  suspended  by  a  wire  that  was  attached  to  the  specimen 
with  an  epoxy  resin  cement,*  the  consistency  and  quick  setting  of  which 
made  it  an  excellent  material.   The  specimen  was  suspended  in  such  a 
manner  that  the  air  stream  always  flowed  parallel  to  the  exposed  cylindri- 
cal surface.  The  top  and  bottom  were  coated  with  the  epoxy;  thus  ft  was 
assumed  moisture  was  being  removed  only  from  surfaces  parallel  to  the 
air  stream. 


*  Concresive  No.  1201,  Adhesive  Engineering  Co.,  1^11  Industrial  Rd., 
San  Carlos,  Calif. 
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The  apparatus  that  was  used  for  vacuum  drying  Is  illustrated  in 
Figure  3>   It  consisted  of  a  glasa  chamber  that  contained  the  sample  and 
that  could  be  evacuated  while  it  was  suspended  from  the  arm  of  an  analyt- 
ical balance. 

Specimen  preparation  was  similar  to  that  used  in  the  air-drying  runs 
except  for  the  attachment  of  the  wire.  The  specimen  was  floated  on  a 
small  pool  of  mercury  in  an  effort  to  provide  a  rapid  ingress  of  heat 
partially  to  compensate  for  the  cooling  of  the  sample  due  to  the  evapora- 
tion of  the  water. 

The  vacuum  vessel  remained  attached  to  the  analytical  balance  through- 
out the  test.  At  four  minute  intervals  the  stopcocks  at  point  A  and  B 
were  sequentially  closed.  The  vacuum  line  was  immediately  disconnected 
at  point  C,   isolating  the  vacuum  vessel  for  the  weight  determination.   The 
time  interruption  averaged  about  45  seconds.  The  vacuum  line  was  then 
without  delay  reconnected  to  the  vessel:   stopcock  B  was  opened  after 
allowing  15  seconds  to  evacuate  the  small  volume  between  A  and  B.  A 
mechanical  pump  was  used  and  the  vacuum  achieved  was  about  0.2  tcrr.  A 
possible  influence  of  the  frequency  of  interruption  for  weighings  was 
checked  and  found  to  be  not  important. 

Mercury  Porosimetry 

Mercury  porosimetry  was  selected  as  the  best  conventional  method  for 
measuring  pore  size  distributions  of  porous  media. 

The  underlying  principle  of  mercury  intrusion  is  theoretically  simple, 
and  revolves  around  the  relationship  of 

r„_22^£SJ»  (6) 
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Fig.    3      Vacuum    Drying    Apparatus 
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where  r  is  the  pore  radius,  P  the  capillary  pressure,  Y   the  surface 
tension,  and  9  the  contact  angle  of  the  three-phase  interface.  Accord- 
ingly, in  mercury  porosimetry  one  determines  the  pore  size  by  calculation 
using  Equation  ( 6  )  and  the  measured  intrusion  pressures  that  are  applied 
incrementally  in  the  apparatus. 

A  source  of  error  in  the  use  of  Equation  ( 6  )  is  the  assumption  that 
the  pores  have  a  circular  cross  section.   If  other  pore  shapes,  e.g. 
rectangular,  ellipitical,  or  triangular,  were  assumed,  the  coefficient 
would  change,  but  this  shift  would  not  change  the  shape  of  the  pore  size 
distribution  curve.  Likewise,  the  value  of  r  depends  upon  the  values 
assumed  for  the  surface  tension  and  contact  angle,  but  once  again,  the 
shape  of  the  distribution  curve  would  not  be  affected.   For  this  study, 
the  values  of  Ritter  and  Drake  (24)  were  used  (  Y  m   480  dynes/cm  and 
9  "  140  ) .  These  values  are  the  conventional  ones  that  have  most 
frequently  been  used  by  others. 

The  Portland  Cement  Association,  a  few  years  before  adequate  commer- 
cial equipment  was  available,  authorized  Landgren  (50)  to  design  and 
build  a  high-pressure  mercury  porosimeter  to  fulfill  the  research  require- 
ments of  the  Association.  The  apparatus  is  pressured  hydraulically  and 
has  a  maximum  working  pressure  of  45000  psl,  which  corresponds  to  a  mini- 
mum pore  that  can  be  intruded  of  about  45  Angstrom  units  in  diameter. 
The  dilatometer  sample  chamber  is  constructed  of  stainless  steel  and  the 
volume  changes  are  determined  by  a  magnetic  sensing  arrangement.  The 
volume  change  accuracy  of  the  system  is  t  0.0001  ml  for  the  usual  size 
of  dilatometer. 

The  sample,  crushed  to  8-16  mesh,  was  placed  in  the  dilatometer  and 
evacuated  with  a  mercury  diffusion  pump  to  10   torr.  The  dilatometer  was 
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then  filled  with  mercury  and  placed  in  the  pressuring  cell.  The  technique 
of  filling  the  dilatometer  was  such  that  the  enclosed  sample  was  subject 
to  an  initial  pressure  of  19-1  psi,  i.e.  pores  larger  than  11 A  microns 
in  diameter  were  already  intruded  at  the  start  of  the  run  and  so  do  not 
appear  in  the  results  of  these  determinations. 

Pressure  was  then  applied  in  increments  to  obtain  equal  length 
changes  in  the  dilatometer  (length  changes  are  directly  proportional  to 
volume  displacements).  The  final  result,  then,  was  a  relationship  between 
applied  pressure  and  volume  of  pore  space  intruded  by  the  mercury.  These 
data  were  then  operated  on  in  the  usual  way  (24)  to  obtain  the  pore  size 
distribution  curve. 

Tortuosity 

Tortuosity,  K.  ,  is  defined  as  the  square  of  the  ratio  of  the  tortuous 
length  (the  average  path  length  a  fluid  "particle"  travels  between  two 
points),  L  ,  to  the  straight  line  length,  L,  between  the  points.  This 
parameter  was  reported  by  Carman  (18)  in  his  treatment  of  the  Rozeny 
equation  and  permeability.  The  measurement  of  tortuosity  poses  difficulty. 
Several  workers  (  51»  52,  53)  have  used  electrical  methods.  Wyllie  and 
Spangler's  (53)  work  was  perhaps  the  most  significant.  Dolch  (8  ) 
measured  tortuosity  on  several  concrete  aggregates  by  electrical  resistance 
methods,  but  expressed  some  reservation  about  the  method. 

The  validity  of  the  electrical  analogue  lies  in  the  assumption  that 
the  electrical  path  in  the  conducting  fluid  that  occupies  the  pores  of 
the  porous  medium  is  the  same  as  the  tortuous  path  of  a  flowing  fluid. 
Furthermore,  It  is  assumed  that  the  solid  material  has  an  infinite  resist- 
ance, and  to  reduce  the  error  of  this  assumption,  a  fluid  with  a  high 
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electrical  conductivity  is  used,  i.e.,  5  percent  sodium  chloride 
solution. 

The  formation  factor,  F,  is  defined  as  the  ratio  of  the  electrical 
resistance  of  the  saturated  porous  material,  R  ,  to  the  electrical  resist- 
ance, R  ,  of  a  body  of  the  saturant  fluid  of  the  same  dimensions. 

F-£  (7) 

Re 

The  formation  factor  is  related  to  the  tortuosity,  K  ,  by 

Kt  -  (Fe)2  (8) 

For  most  porous  materials,  the  pore  space  Is  only  partially  saturated 
with  a  liquid  (even  under  vacuum  saturation  conditions),  causing  the 
formation  factor,  F  ,  to  be  high.  The  following  relationship  then  applies 

(53) 

i  -  s  "n  (9) 

F       w 

where  S  is  the  degree  of  saturation,  and  the  exponent,  n,  has  been 
shown  (52)  to  be  equal  to  2  -  1/2  until  a  relatively  low  saturation  pre- 
vails. 

Cylindrical  samples  were  vacuum  saturated  with  a  5  percent  (50  g 
in  1000  ml  of  solution)  NaCl  solution.  After  2^  hours  of  immersion  they 
were  surface-dried  on  a  towel  slightly  wet  with  the  solution.  The  ends 
were  then  covered  with  graphite  paste,  and  the  resistance  was  measured 
using  a  1000  cps  alternating  current  bridge.  The  specific  resistance  of 
the  solution  was  measured  independently,  and  the  result  agreed  well  with 
values  reported  in  standard  handbooks.  The  values  of  tortuosity  were 
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then  calculated  from  the  data  by  the  use  of  Equations  (8  and  9  )  and  the 
assumption  that  n  ■  2. 

Tensile  Strength 

The  tensile  splitting  test  was  selected  as  the  best  procedure  for 
measuring  the  tensile  strength  of  the  aggregate.  Three  cylindrical  cores, 
1  in.  long  and  0.67  in.  in  diameter  were  carefully  selected  to  eliminate 
any  visible  fractures  in  the  specimens  and  were  oven  dried.  The  aggregate 
designated  IDE  was  so  highly  fractured  naturally  that  no  specimen  was 
found  to  be  sound  enough  for  such  testing. 

A  supplementary  loading  plate  had  a  ball  and  joint  assemblage  tbat 
allowed  free  movement  for  alignment.  Two  bearing  strips  of  cardboard 
approximately  1/16  in.  thick  and  1/2  in.  wide  and  having  a  length  slightly 
longer  than  the  specimen  were  placed  between  the  specimen  and  the  bearing 
plates.  The  loading  rate  was  600  lb/min.  using  a  hydraulic  universal 
testing  machine. 

Freezing  and  Thawing  Tests 

Concrete  prisms,  3  x  3  x  11  in* >   were  prepared  and  were  tested  in 
accordance  with  ASTM  Method:  C-291  (5*0. 

Ledge  rock  was  crushed  to  3A-in.  maximum  size  with  a  laboratory 
crusher.  These  materials,  the  gradations  of  which  are  listed  in  Table  3, 
were  then  vacuum  saturated  with  water  and  brought  to  a  saturated  surface  - 
dry  state  just  prior  to  mixing.   Sand  from  Elgin,  Illinois  was  used  as 
the  fine  aggregate.   Its  gradation  is  given  in  Table  h. 

The  concrete  was  designed  to  provide  optimum  conditions  that  would 
delineate  the  behavior  of  the  coarse  aggregate.  The  properties  of  the 
mixes  are  given  in  Table  5«  The  mixes  were  designed  to  have  a  constant 


Table  3 
Coarse  Aggregate  Gradation 


Sieve 

Percent 

Retained 

Percent 

No. 

Individual 

Cumulative 

Passing 

3/4« 

0 

0 

100 

3/8" 

63 

63 

37 

4 

37 

100 

0 
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Table  4 

Fine 

Aggregate 

Gradation 

Sieve 

Percent  Retained 

Percent 

No. 

Individual 

Cumulative 

Passing 

4 

0 

0 

100 

8 

18 

18 

82 

16 

15 

33 

67 

30 

24 

57 

43 

50 

30 

87 

13 

100 

8 

95 

5 

Pan 

5 

100 

0 

Table  5 
Concrete  Mix  Properties 
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Id.  No. 

Aggregate 
Cement  Ratio 
Bv  Weight** 

3.12 

Unit  Weight      Air* 
Lb/Ft. 3       Percent 

148.9        4.6 

Slump* 
Inches 

1-1S 

3 

6?-2S 

Material  Not  Available 

4?-2S 

2.92 

5.6 

3 

9-1S 

2.86 

4.8 

3 

MM 

Material  Not  Available 

ML 

2.99 

148.2         5.^ 

1 

MTS 

3.00 

146.0         5.4 

2 

IS 

3.01 

145.0        5.3 

3  1/2 

II 

2.86 

143.5        5.9 

3 

IDO 

2.94 

145.3        6.0 

3 

IDE 

2.93 

143.4        6.2 

3 

LA 

2.90 

141.2        7.4 

3 

*  Due  to  a  limited  supply  of  aggregate,  batch  sizes  were  limited  to 
0.150  cubic  feet.  Air  content  was  measured  with  a  500  oc.  mortar 
pressure  meter.  Slumps  were  estimated. 

**  Aggregate  volumes  were  kept  constant. 
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cement  factor  (6  bags  per  cubic  yard),  a  constant  aggregate  content  by 
volume  (66.7  percent),  a  uniform  workability,  and  an  air  content  that 
would  ensure  freedom  of  the  paste  from  damage  under  freezing.  Therefore 
any  frost  damage  in  these  specimens  was  assumed  to  be  caused  by  the 
aggregates  under  test. 

Two  aggregates,  67-2S  and  MM,  were  not  included  in  this  portion  of 
the  study  because  of  insufficient  material.  However,  several  investigators 
(  21»  55    )  have  reported  work  with  67-2S,  and  there  is  general  agreement 
that  its  durability  record  is  excellent. 

After  lU  days  cure  immersed  in  saturated  lime  water,  the  freezing 
and  thawing  exposure  was  begun.  Approximately  6  cycles  per  day  were 
experienced.  The  freezing  was  in  air  at  0  F,  and  the  thawing  was  in  water 
at  kO   F. 

From  time  to  time  measurements  were  made  of  length  change  to  the 
nearest  0.0001  in.,  of  weight  change,  and  of  sonic  modulus  of  elasticity, 
using  conventional  equipment. 
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Data 


General 

Average  values  for  bulk  density,  porosity,  absorption,  degree  of 
saturation,  tensile  strength,  and  tortuosity  are  reported  in  Table  6. 
The  volume  of  voids,  used  to  caloulate  bulk  density,  porosity,  and  degree 
of  saturation  was  measured  with  a  McLeod  porosimeter  (21).  The  degree  of 
saturation  is  the  ratio  between  the  void  volume  filled  with  the  saturant 
and  the  total  void  volume.  The  bulk  volume  of  the  solid  was  determined 
from  the  external  dimensions  of  the  specimen.  Tortuosity  values  were 
calculated  by  equations  8  and  9» 

The  aggregates'  insoluble  residue  content  and  their  qualitative 
identification  by  X-ray  diffraction  are  given  in  Table  7.  A  ten  percent 
solution  of  hydrochloric  acid  was  used  to  dissolve  the  carbonation 
components  from  the  aggregates. 

Evaporation  Rates 
The  evaporation  data  for  vacuum  and  air  drying  were  collected 
initially  with  weight  change  as  a  function  of  time.  The  data  for  vacuum 
drying  are  tabulated  in  Appendix  D.  Equation  G9  suggests  that  the  rate 
of  evaporation  during  the  second  falling  rate  period  is  inversely 
proportional  to  the  accumulated  volume"  of  fluid  loss.  Shown  in  Figure  4 
are  a  few  typical  drying  rate  curves.  Observe,  that  the  first  few  points 
lie  in  the  first  falling  rate  period,  and  are  not  needed  to  determine  the 
slope  of  the  curve  in  the  second  falling  rate  period.  Thus,  for 
simplicity,  the  data  that  lie  only  in  the  second  falling  rate  period 
for  the  remaining  specimens  are  shown  in  Figs.  5  to  8. 
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Table  6 

Density,  Porosity,  Absorption,  Tensile  Strength, 
Saturation  and  Tortuosity 


Id.   No. 

Bulk 

Porosity 

Absorption 

Saturation** 

Tortuosity** 

Tensile 

Density, 
^b-g/cra^ 

€ 

* 

Sw 

Kt 

Strength 

pt 

PS  I 

1-1S 

2.80 

0.029 

0.7 

0.89 

1000 

2200 

67-2S 

2.60 

o.d+o 

1.1 

0.69 

56 

1600 

1*7- 2S 

2.  1*2 

0.105 

k.6 

O.89 

60 

1550 

9- is 

2.28 

0.206 

9.8 

1.00 

8.1 

700 

ML 

2.35 

0.117 

5-0 

0.98 

& 

800 

MTS 

2.67 

0.01*7 

1.8 

O.96 

196 

950 

MM 

2.38 

0.105 

^•5 

0.62 

23 

11*00 

IS 

2.79 

0.056 

1-5 

0.79 

1*68 

1600 

II 

2.1*6 

0.079 

3-2 

0.92 

118 

1150 

IDE 

2.61* 

0.031* 

1.7 

0.83 

325 

— 

LDU 

2.63 

0.075 

2.0 

0.73 

75 

1150 

IA 

2.58 

0.051* 

1-5 

0.80 

82 

11*00 

FGF 

lAl 

0.301* 

21.6 

1.00* 

__ 

__ 

FGM 

I.58 

0.256 

16.0 

1.00* 

__ 

— 

FGC 

1.1+6 

0.218 

ll*. 9 

1.00* 

-- 

— 

FB 

5.1*6 

0.300 

l*.5 

1.00* 

« 

-- 

Vycor 

1.20 

0.357 

21.7 

0.86 

10.6*** 

-- 

w/c-.35 

1.86 

0.281* 

15.2 

1.00* 

10.9 

-- 

w/c».65 

1.2l* 

0.1*65 

37-5 

-- 

-- 

-- 

w/c-.8o 

1.19 

0.602 

1*6.6 

-- 

5-9 

-- 

Brick 

0.163 

7.0 

*  These  values  were  assumed. 
**  These  values  were  measured  on  different  specimens. 
***  See  Appendix  A. 
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Table  7 

Insoluble  Residue  Contents  and  Qualitative 
Identification  by  X-ray  Diffraction 

Insoluble  Residues 
Id.  No. 

Percent       Composition 


1-1S 

3.38 

Quartz,  mica 

67.2s 

2.90 

Quartz,  mica 

47-2S 

9.07 

Quartz,  mica 
Kaolinite 

9-1S 

7.19 

Quartz,  mica 

ML 

1.71 

Quartz,  mica 

MTS 

1.80 

Quartz 

IS 

2.01 

Quartz,  mica 

II 

0.5^ 

Quartz,  some  clay 

I  DO 

1.11 

Mica,  clay 

I  DE 

5.52 

Quartz,  mica 
Montmorillonite 

I  A      0.59      Kaolinite,  mica,  chlorite 
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In  Table  8  are  given  the  porosity,  exposed  area  to  evaporation, 
tortuosity,  and  the  slope  of  the  curves  found  in  Figs.  4  to  8  for  those 
aggregates  that  were  dried  by  the  vacuum  technique.  These  properties 
were  all  measured  on  the  same  specimen  for  each  aggregate  material. 

The  data  for  those  specimens  that  were  air-dried  are  shown  in  two 
sets  of  curves.  Figs.  9  to  15  show  the  logarithm  of  the  drying  rate  as 
a  function  of  moisture  content.  These  curves  are  of  interest  because  of 
their  similarity  to  those  presented  by  other  investigators  (39»  ^0« 

The  above  rate-mass  curves  have  a  parabolic  shape  in  the  first 
falling  rate  period  for  about  the  first  40  percent  of  the  moisture  loss  — 
the  exact  value  depends  upon  the  material's  pore  structure  and  surface 
environment.  If  these  data  are  presented  as  the  logarithm  of  accumulated 
volume  less  as  a  function  of  the  logarithm  of  time,  a  straight  line 
relationship  will  exist  for  that  portion  of  the  previous  curve  that  was 
parabolic.  The  mathematical  relationships  between  these  two  functions 
are  shown  in  Appendix  B. 

Figs.  16  to  18  show  the  data  of  the  air-dried  specimens  in  a  form 
that  is  similar  to  the  data  on  vacuum  drying  presented  in  Fig.  4. 

Mercury  Porosimetry 

The  pore  size  distribution  curves  determined  by  mercury  intrusion 
are  shown  in  Figs.  19  to  25.  The  median  pore  size,  calculated  from  the 
logarithmic  distribution  is  shown  on  the  respective  figure. 

The  distribution  function,  the  ordinate  of  the  distribution  curve, 

is  the  differential  of  the  accumulative  volume-size  function.  It  has  the 

2   3 
dimension  of  cm  /cm  ,  and  is  proportional  to  the  internal  surface  area 
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Table  8 

Slope,  "b",  for  Vacuum  Dried  Aggregates 
and  Their  Respective  Forosity,  Tortuosity  and  Area 


Id  No. 

Porosity 

Area2 
A,   cm 

Tortuosity 
Kt 

Slope              0 
b,  gr-cm/min  x  10 

1-1S 

- 

- 

- 

6?-2S 

3.6 

11.00 

4.0 

0.118 

47- 2S 

10.8 

13.15 

7.3 

0.219 

47-2S(o) 

10.1 

10.10 

7.7 

0.288 

9-1S 

21.3 

14.21 

3.3 

1.62 

MM(1) 

11.5 

13.40 

7.4 

1.559 

MM(2) 

15.1 

13.70 

5.1 

2.613 

ML(1) 

14.6 

12.81 

4.7 

O.833 

ML(2) 

13.6 

12.82 

6.4 

0.762 

ML(3) 

15.2 

12.73 

5.6 

1.29 

ML(4) 

13.O 

13.05 

6.5 

1.39 

MTS(l) 

6.6 

13.00 

12.6 

0.386 

KTS(2) 

6.2 

12.65 

16.0 

0.328 

IS 

5.2 

14.00 

23.8 

0.014 

IKD 

10.5 

13.50 

4.8 

0.765 

11(2) 

6.8 

13.18 

11.4 

O.438 

IDU(l) 

6.7 

13.88 

14.7 

0.497 

IDU(2) 

12.2 

13.00 

7.8 

2.585 

IDE 

4.8 

12.85 

10.3 

0.097 

IA(1) 

5.1 

13.50 

8.3 

0.079 

IA(2) 

6.4 

13.13 

6.5 

0.094 
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at  "d".  The  area  under  the  curve  between  pore  sizes,  cL  and  &2*   ^s 
equivalent  to  the  volume  of  those  pores. 

Freezing  and  Thawing 
The  concrete  prisms,  tested  in  accordance  to  ASTM  Designation  C-291, 
"Resistance  of  Concrete  Specimens  to  Rapid  Freezing  in  Air  and  Thawing  in 
Water,"  were  measured  for  length  change,  weight  change,  and  relative 
dynamic  modulus  of  elasticity.  The  percent  expansion  data  are  shown  in 
Figs.  26  to  29,  while  the  percent  weight  change,  relative  dynamic 
modulus  of  elasticity,  and  photographs  of  the  tested  specimens  are  shown 
in  Appendix  E. 
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DISCUSSION 

Theoretical  Considerations 
This  discussion  will  be  limited  primarily  to  those  concepts  that 
apply  directly  to  evaporating  phenomenon.   The  drying  process,  which  was 
treated  in  some  detail  in  previous  sections  of  this  text,  can  be  character- 
ized generally  by  three  regimes:   l)  the  constant  rate  period,  2)  the 
first  falling  rate  period  and  3)  the  second  falling  rate  period.  A  brief 
review  of  each  of  these  regimes  will  permit  the  selection  of  the  drying 
regime  that  most  accurately  characterizes  the  pore  systems  of  consolidated 
porous  media. 

Constant -Rate  Period 
Drying  during  the  constant-rate  period  is  a  process  whereby  the  porous 
medium  is  able  to  maintain  an  evaporating  surface  that  is  equivalent  to  a 
surface  of  free  water.   For  this  state  to  prevail,  an  adequate  source  of 
fluid  must  exist  and  it  must  be  able  to  present  itself  at  the  surface  at 
a  rate  that  is  at  least  commensurate  to  the  rate  of  evaporation  from  a 
free  water  surface.   The  data  show  that  this  requirement  is  not  met  by 
most  of  the  samples  studied  here.  The  matter  is  discussed  in  more  detail 
later. 

First  Falling  Rate  Period 
When  the  exposed  surface  is  unable  to  receive  sufficient  fluid  to 
supply  the  environmental  demand,  the  first  falling  rate  period  commences; 
the  surface  begins  to  dry,  and  its  effective  area  of  evaporation  decreases. 


6k 

This  Is  the  region  in  which  the  liquid  stays  preferentially  in  the  smaller 
pores  because  of  their  higher  capillary  potential  and  in  which  the  larger 
pores  become  increasingly  emptied. 

As  the  fluid  recedes  in  the  larger  pores,  the  diffusion  process 
begins  to  influence  the  rate  of  evaporation.  The  first  falling  rate 
period  is  a  transition  zone  (kk,   kf)    in  which  the  fluid  lost  from  the 
surface  of  the  sample  becomes  progressively  less,  and  the  fluid  lost  by 
diffusion  from  the  interior  becomes  progressively  greater. 

Second  Falling  Rate  Period 

Eventually  the  drying  becomes  predominately  a  diffusional  process,  a 
sub-surface  evaporation  in  which  a  quasi-plane  of  vaporization  retreats 
farther  and  farther  into  the  porous  medium. 

By  considering  a  long  circular  cylinder  in  which  diffusion  is  radial, 
and  the  concentration  is  a  function  of  the  radius,  r,  and  time,  t,  only, 
and  by  employing  cylindrical  coordinates,  equation  (2)  becomes 

Following  essentially  the  treatment  by  Crank  (^6), 

C  -  -i>  exp  (-DoC2t)  (11) 

is  a  solution  for  equation  (10)  provided  D  is  constant  and  ^  is  a  function 
of  r  that  will  satisfy  Bessel's  equation  of  order  zero,  which  is 

£*+k    ±3>+flt8V-0  (12) 

.  2   r  dr  %   ' 

dr 

Applying  the  appropriate  boundary  conditions (C  ■  f  (r),  0  <  r(i, 
t  «  0;  C  ■  0,  r  «  a,  t  >  0)  Crank  (k6)   gives  as  the  solution 


o© 


k 


"~    Hl^T^    e*P(-D°C2nt)      (13) 


n"! 
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where  W  is  the  average  moisture  content  at  time  t,  and  the  OC    's  are 
positive  roots  of 

J  (a  ♦  OC    )  -  o  (Ik) 

o       n'  v   ' 

and  J  (x)  is  the  Bessel  function  of  the  first  kind  of  order  zero.   Some 

ox  ' 

of  the  contingent  restrictions  and  assumptions  inherent  in  equation  (13) 
are  discussed  in  Appendix  I. 

During  the  later  stages  of  drying  by  vacuum  it  is  assumed  that  the 
drying  rate  is  controlled  by  flow  of  vapor  through  the  porous  medium  and 
that  this  movement  is  described  by  the  equation  of  slip  flow. 

In  Appendix  G  is  derived  equation  (G9),  which  describes  the  flow 
process  of  the  vapor  during  the  drying  of  a  wet  porous  solid.   It  will 
be  observed  that  this  equation  is  linear  in  the  rate  of  weight  change  of 
the  sample  and  reciprocal  volume  of  saturant  lost.  The  slope  of  the 
curve  contains  only  the  hydraulic  radius  of  the  pores  and  quantities  that 
are  known  or  can  be  reasonably  assumed. 

One  of  the  assumptions  in  the  derivation  of  equation  (G9)  is  that 
the  liquid  recedes  into  the  pores  as  the  drying  process  proceeds  and  in 
so  doing  leaves  no  residue  on  the  walls,  I.e.  that  at  any  time  the  porous 
body  can  be  divided  into  a  still-completely  saturated  region  and  one 
completely  dry  the  pores  of  which  contain  only  the  flowing  vapor. 

This  assumption  is  only  very  roughly  true,  as  was  shown  by  experi- 
ments in  which  the  drying  process  was  interrupted,  the  sample  cleaved 
into  several  fragments,  and  the  moisture  content  of  these  determined.   Since 
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the  walls  do  remain  wet  during  the  drying  process  the  actual  average 
flow  path  is  shorter  than  that  envisaged  in  the  aforementioned  assumption 
and  embodied  in  equation  (G9)  and  it  is  then  necessary  to  correct  L  in 
equation  (G7)  to  a  smaller  value.  Intuitively  one  would  expect  this 
correction  to  be  a  function  of  some  physical  property  of  the  porous 
medium.  The  degree  of  wall  wetness  and  thus  the  effective  reduction  of 
the  flow  length  are  probably  functions  of  specific  surface  area,  S,  and 
pore  size.   Ihe  porosity  is  related  to  these  parameters  by 

e  =  Sra  (15) 

Thus  if  it  is  assumed  that  the  effective  flow  length 

Leff.  -  eL  (l6) 

and  this  is  incorporated  in  equation  (G7)  the  result  is 

2.       2—   h  S.m 
dW   eA>!AP 


tl±£  +  1 )  I  (17) 


dt   Kt1/2RT   V?     *1     V 

by  means  of  which  the  value  of  m  can  be  calculated  from  the  data. 

Permeability 
The  coefficient  of  permeability  which  is  usually  referred  to  as 
simply  "permeability"  expresses  the  resistance  to  flow  experienced  by  a 
fluid  in  a  porous  medium  and  is  defined  by 

2-  .  Ijf  (18) 

At    AL 

where     V  =  volume  of  fluid  flowing  through  the  medium  in  time  t 

A  =  approach  or  macro  area  of  medium 

K  =  permeability 

1  as  viscosity  of  fluid 
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AP  *  head,  in  pressure  units,  between  inlet  and  outlet  faces 

AL  ■  overt  flow  path 

2 
The  cgs  unit  of  permeability  is  the  cm  .  The  conversion  into  other 

2  11 

commonly-used  units  is  1  cm  =  1.013  x  10   millidarcy  (md)  ■  1.097  x 

10  cm/sec  (if  the  fluid  is  water  at  room  temperature). 

Results 

The  following  results  have  been  calculated  using  the  concepts  and 
relationships  developed  in  the  previous  section. 

The  average  pore  sizes  calculated  by  three  methods  are  shown  in 
Table  9«  One  pore  diameter  was  calculated  from  equation  (17)  using 
"vacuum  dried"  data.  The  data  used  are  on  Figs,  k   to  8.  The  pore 
diameters  calculated  from  mercury  porosimetry  are  the  logarithmic  means 
of  the  curves  shown  in  Figs.  19  to  25-   It  is  recognized  that  the  choice 
of  a  logarithmic  mean  is  arbitrary.  Had  the  algebraic  mean  been  used  the 
resultant  pore  size  would  be  unreal istically  large.  A  further  basis  for 
this  choice  is  the  reasonable  permeability  calculated  from  the  logarithmic 
increment  method  (q.v.).  The  data  for  the  direct  permeability  measure- 
ments is  from  Dolch  (8). 

The  pore  size  of  a  porous  medium  can  be  related  to  its  permeability 
by  the  Kozeny-Carman  equation  (l8), 

K-#-  (19) 

O  t 

If  the  pore  diameters  calculated  from  the  evaporation  data  are  used 
in  this  equation,  with  the  appropriate  porosities  and  tortuosities  from 
Table  8,  permeabilities  are  obtained  and  reported  in  Table  10, 
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Iden. 
No. 

1-1S 
6T-2S 
1+7-2S 

9-1S 


MM 

(1) 

MM 

(2) 

ML 

(1) 

ML 

(2) 

ML 

(3) 

ML 

(4) 

MTS 

(i) 

MTS 

(2) 

IS 

II 

(1) 

II 

(2) 

IDU 

(1) 

IDU 

(2) 

IDE 

IA 

(1) 

IA 

(2) 

Table  9 

Calculated  Pore  Diameter 

Pore  Diameter,  microns 
Mercury 
Evaporation    Porosimetry    Permeability  (  8 ) 


0.35 
0-97 
2.93 
6.23 
6.90 
2.38 
2.93 
4.27 
5-11 
2.88 
3.18 
0.19 
2.13 
2.92 
3.98 
10.41 
O.627 
0.39 
O.38 


0.77 
0.08 

0.35 
I.85 

0.46 

0.59 


0.71 
0.10 

0.37 

0.40 

0.32 
0.31 


0.64 
0.52 
0.46 
I.58 
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Table  10 

Calculated  Permeabilities 

Iden. 
No. 

Evaporation 

Permeability,  md 
Mercury 
Porosimetry    Permeability  (  8 ) 

1-1S 

0.0012 

0.00205 

67-2S 

0.0064 

0.00014 

0.0072 

47-2S 

0.039 

0.023 

0.032 

9-ls 

2.0 

1.48 

4.30 

MM  (1) 

1.89 

0.18 

MM  (2) 

6.4o 

ML  (1) 

0.872 

0.12 

ML  (2) 

0.695 

ML  (3) 

2.04 

ML  (4) 

1.87 

MTS  (1) 

0.798 

0.023 

MTS  (2) 

O.568 

IS 

0.00073 

0.0006 

II  (1) 

0.48 

0.005 

II  (2) 

0.10 

IDU  (1) 

0.12 

0.59 

n>u  (2) 

5.03 

IDE 

0.043 

0.012 

U  (1) 

0.026 

0.024 

IA  (2) 

0.051 
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assuming  a  shape  factor,  K  ,  of  2.7  (53).  The  permeabilities  of  the 
Indiana  aggregates  are  also  reported  (8). 

A  slight  modification  of  equation  (19)  will  permit  permeabilities  to 
be  calculated  from  the  pore  size  distribution  curves  derived  from  mercury 
intrusion.   The  distribution  curve  is  divided  into  n  sections,  allowing 
each  logarithmic  segment  of  the  abscissa,  log  Am,  to  be  equal  in  size, 
and  then  the  fractional  contribution  of  each  section  to  the  total  permea- 
bility is  computed.   The  permeability  for  one  such  section  can  be  written 
as 

h-fr    *i*a«i2  (20) 

o  t 

where  the  product  of  Y.  A  m^  is  the  fraction  of  porosity  between  m.  and 

m?.  The  total  permeability  is  the  summation  of  those  of  the  individual 

segments 

n 

K  "  fV  £  Yi  ^  Am  (21) 

oKt  \%1  i 

Using  equation  (21)  permeabilities  have  been  calculated  from  the  porosimetry 

data  and  are  shown  in  Table  10.   Purcell  (56)  used  a  similar  treatment. 

The  coefficients  of  diffusion  are  given  in  Table  11.  These  results 
were  calculated  from  equation  (13)  and  from  the  slope  of  the  second 
falling  rate  period  in  Figs.  16  to  18.  The  rate,  dW/dt,  and  W  were 
evaluated  using  Figs.  9  to  12. 

It  is  customary  to  calculate  a  durability  factor  for  freeze-thaw 
exposure  from  the  relative  dynamic  modulus  of  elasticity  data  (Figs.  S6 
to  E8).  The  durability  factor,  DF,  is  defined  in  ASTM  Method:  C-291  (5*0 
and  is  a  relative  measure  of  the  durability  of  the  concrete  in  this  test, 
the  higher  the  better.  The  calculated  durability  factors  are  given  in 
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Ident. 

No. 


Table  11 
Coefficients  of  Diffusion 


Coefficient  of  Diffusion 


cm  /sec  x  10 


1-lS 

0 

,40 

67-  2S 

0.54 

47- 2S 

2.18 

9-1S 

3 

.21 

MM 

' 

1 

.17 

ML 

2 

.10 

MTS 

9 

.16 

IS 

0 

.48 

IA 

0 

.7? 

II 

2. 

.22 

IDC 

1 

.17 

IDE 

1, 

.62 

Table  12 

Freezing  and 

Thawing  Results 

Ident 

• 

Durability 

Expansion 

Failure 

No. 

- 

Factor 

Factor 

Type 

1-lS 
1-lS 

(A) 
(B) 

75 
62 

8.3 
10.0 

Sound 

67-2S* 

98 

- 

47- 2S 

(A) 

62 

9.1 

Sound 

47-2S 

(B) 

75 

11.0 

9-1S 

(A) 

23 

5.0 

9-1S 

(B) 

22 

5.0 

Boundary 

9-1S 

(C) 

8 

3.0 

ML 

ML 

(A) 
(B) 

17 
17 

1.4 
1.4 

Boundary 

MTS 

(A) 

37 

6.7 

MTS 

(B) 

43 

7.3 

Boundary 

IS 

(A) 

30 

5-3 

Internal 

IS 

(B) 

30 

5.6 

Fracture 

IA 
IA 

(A) 
(B) 

91 
89 

29.4 
37.1 

Sound 

II 

(A) 

65 

9.4 

Sound 

II 

(B) 

54 

11.0 

IDU 

(A) 

81 

10.0 

Sound 

IDU 

(B) 

68 

12.3 

IDE 
IDE 

(A) 
(B) 

59 
64 

8.0 
8.2 

Sound 

*  This  material  was  not  available  for  testing. 
Value  reported  is  from  McLaughlin  {55)» 
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Table  12  together  with  a  similar  expression  for  the  expansion  performance. 
This  expression  is  called  the  Expansion  Factor,  EF,  and  is  defined  as 


EF-^  (22) 


where 

e  =  percent  expansion  at  N  cycles 

i  ■  number  of  cycles  at  which  e  reaches  one-tenth  of  one 
percent  or  300  cycles,  whichever  is  first,  and 

h  ■  specified  number  of  cycles,  300,  at  which  the  exposure 
is  to  be  terminated 

Discussion  of  Data  and  Results 


The  objective  of  this  section  will  be  limited  primarily  to  evaluating 
the  evaporation  test  as  a  technique  for  measuring  the  pore  characteristics 
of  porous  media.  The  section  to  follow  will  interpret  the  results  and  the 
application  of  the  measured  pore  parameters  and  their  relevance  to  the 
durability  of  concrete. 

Evaporation  by  Air  Drying 
Figs.  9  to  13  present  the  drying  data  in  a  form  generally  employed 
in  drying  technology.  From  these  curves,  it  is  apparent  that  there  is  no 
strictly  constant  rate  period  for  the  materials  tested  in  this  investiga- 
tion, although  several  have  a  quasi-constant  rate  region.  For  porous 
materials  that  are  similar  to  aggregates  and  concrete,  other  investigators 
(5?»  58)  have  found  no  constant  rate  period.  A  constant  rate  period  did 
exist  for  the  fritted  glass  specimens,  Fig.  1^.  These  were,  however,  of 
a  different  shape  than  that  of  the  aggregate  specimens;  the  fritted  glass 
samples  were  cut  from  filtering  funnels  and  the  evaporation  was  from  the 
flat  parallel  faces  so  the  flow  path  was  short,  about  1  mm.  This  geometry 
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resulted  in  such  short  flow  paths  that,  it  is  surmised,  the  flow  to  the 
surface  was  easily  enough  to  maintain  a  fairly  constant  evaporation  area 
and,  thereby,  a  constant  evaporation  rate.  This  was  not  true  of  the  other 
samples  where  the  flow  path  was  longer. 

It  can  be  seen  from  the  curves  for  the  performance  of  the  fritted 
glasses  that  the  pore  size  does  not  influence  the  magnitude  of  the  constant 
rate.  The  pore  sizes  varied  by  about  an  order  of  magnitude,  yet  the 
evaporation  rates  were  all  about  the  same. 

For  most  of  the  samples,  except  the  fritted  glasses,  a  first  falling- 
rate  period  occurred.  As  discussed  in  Appendix  B  ,  a  mathematical 
expression  can  be  written  for  this  region  of  drying,  which  describes  the 
parabolic  shape  of  the  curves.  An  analysis  of  the  intercepts  and  slope 
parameters  of  the  curves  plotted  according  to  equation  Bl  shows  no 
correlation  whatsoever  with  any  pore  characteristic  measured  by  mercury 
intrusion  or  permeability.  From  this  apparent  lack  of  correlation  and 
the  fact  that  this  region  is  a  transition  zone  between  the  constant  rate 
period  (if  one  exists)  and  a  period  of  diffusion,  it  becomes  evident  that 
the  pore  size  (and  aggregate  durability,  presuming  a  relationship)  cannot 
be  resolved  by  the  parameters  of  the  first  falling  rate  period. 

Equation  (17)  suggests  that  the  drying  rate  and  the  reciprocal 
accumulated  volume  loss  are  related  linearly.  The  air-dry  data  plotted 
in  this  form  are  shown  in  Figs,  k     and  8  .  It  can  be  seen  from  these 
figures  that  the  first  and  second  falling-rate  periods  are  Indeed  linear 
in  these  variables. 

However,  some  question  does  arise. as  to  the  validity  of  equation  (l?) 
when  applied  to  air-dried  data.  As  pointed  out  in  the  discussion  of 
theory,  air  drying  causes  a  diffusional  flow  of  the  vapor,  the  pressure 
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difference  being  approximately  equal  to  the  vapor  pressure  of  the  liquid. 
Assuming  the  total  pressure  change  to  be  20  mm  Hg,  the  mean  pressure 
would  be  10  mm  Hg.  Making  the  other  appropriate  substitutions  into 
equation  (1?),   the  pore  sizes  calculated  are  one  to  two  order  of  magnitude 
smaller  than  those  measured  by  permeability  or  mercury  intrusion,  because 
basically  the  equation  does  not  apply  to  the  process  in  question.   However, 
if  the  mean  pressure,  P,  is  increased  to  one  atmosphere,  for  which 
justification  is  difficult,  good  agreement  is  had  between  the  pore  size 
calculated  therefrom  and  other  methods.  There  was  no  direct  method  found 
in  this  study  to  relate  air-dry  data  to  a  mathematical  model  and  calculate 
pore  sizes.  An  empirical  relationship,  however,  exists,  as  shown  above. 
The  data  and  curves  are  included  primarily  to  assist  others  who  may  develop 
alternate  interpretations. 

The  air-drying  rates  in  the  second  period  are  related  to  a  diffusion 
coefficient.  Equation  (13)  was  used  to  calculate  these  coefficients,  and 
the  results  are  given  in  Table  11  .  While  there  is  a  rough  correlation 
between  these  coefficients  and  the  calculated  permeabilities,  it  is  not  a 
good  one.  The  diffusion  coefficient  tends  to  hug  a  central  point,  and 
does  not  differentiate  the  pores  at  either  the  high  or  small  size  range. 

Evaporation  by  Vacuum  Drying 
During  the  second  falling  rate  period  of  vacuum  drying  the  fluid  is 
being  removed  by  a  process  that  is  considered  to  be  described  by  equation 
(  1?  ) .  The  expression  for  the  slope  of  the  curve  contains   l)  parameters 
that  are  measured  experimentally:   e,  A,  K  ,  2)  those  that  are  known  with 
certainty:  M  ■  18  for  water  vapor,  R,  the  gas  constant,  3)  those  that 
can  be  reasonably  assumed:  AP  ■  vapor  pressure  of  water  at  temperature  T, 
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P  ■  AP/2,  17  ■  viscosity  of  water  vapor  at  I,  v  ■  the  mean  molecular  speed 
at  T,  and  T  itself,  which  was  assumed  to  be  room  temperature  (295  K)> 
although  it  is  recognized  that  the  evaporative  cooling  probably  lowered 
the  actual  temperature  in  spite  of  the  thermostatic  effects  of  the  mercury 
pool  on  which  the  sample  floated,  and  k)    those  that  could  be  assumed,  but 
with  perhaps  less  certainty:   K  -  2.7  from  the  work  of  Wyllie  and 
Spangler  (53)  >  K,  =  1  from  Carman  (18),  and  8.  =  0.75>  also  from  Carman. 

The  measured  slopes  were  then  used  to  calculate  the  values  of  the 
hydraulic  radii,  m.   The  results  are  shown  in  Table  9  along  with  those 
determined  by  the  mercury  porosimetry  experiments  and  by  direct  measure- 
ment of  permeability  by  Dolch  ( 8  )  on  the  same  rocks. 

A  further  point  must  be  considered  with  respect  to  equation  (17). 
The  equation,  as  it  stands,  is  of  a  curve  that  passes  through  the  origin, 
which  the  curves  of  the  data  obviously  fail  to  do.  All  have  an  intercept 
on  the  reciprocal  volume  axis.  The  explanation  is  that  V  is  a  limited 
variable,  i.e.  there  is  a  maximum  volume  lost,  which  is  equal  to  the  pore 
volume.  So  there  is  a  minimum  1/V  at  which  the  rate  should  be  zero.  This 
is  the  case  for  all  the  curves.  The  situation  amounts  to  an  axis  trans- 
formation so  that  the  V  of  equation  (17)  becomes 


1  m   1 1 

V  =  V  ,,  "  V 

ef f    max 


(23) 


and  equation  (l?)  becomes 


*  -constant  (^—  -  1— )  (2k) 

eff    max 

It  was  observed  that  the  water  content  of  the  samples  at  the  time 

the  evaporation  rate  dropped  effectively  to  zero,  was  about  a  saturation 

of  5-15  percent.  This  means  that  equation  (17)  fails  to  describe  the 
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very  dry  end  of  the  process,  but  the  pore  size  calculated  from  the  slopes 
of  equation  (17)  should  be  accurate  if  that  equation  is  valid  in  the 
region  over  which  the  data  were  determined. 

Agreement  between  the  pore  diameters  calculated  from  vacuum  evapora- 
tion data  and  those  calculated  from  direct  measurement  of  permeability  is 
excellent.  A  similar  agreement  is  had  between  the  average  pore  diameters 
calculated  from  mercury  intrusion  and  from  evaporation  data  for  materials 
1-1S,  47-2S,  9-1S,  IS,  IDE  and  IA.  For  the  remaining  aggregate  samples, 
the  pore  sizes  measured  by  mercury  porosimetry  were  much  smaller  than 
those  measured  by  evaporation.  This  difference  is  not  alarming  and 
should  be  expected,  as  mercury  porosimetry  measures  the  size  of  the 
entrance  to  the  pore,  and  if  the  porous  medium  has  a  system  of  so-called 
"ink-bottle"  pores,  mercury  porosimetry  underestimates  their  sizes.  Thus 
67-2S,  MM,  ML,  MTS,  II,  and  IDU  may  have  an  "ink  bottle"  pore  system. 

The  permeability  is  readily  calculated  from  pore  6ize  by  means  of 
equation  (19) •  The  permeabilities  so  calculated  from  pore  sizes  determined 
by  evaporation  and  mercury  porosimetry  are  shown  in  Table  10  .  Agreement 
between  the  two  is  good,  including  also  a  comparison  with  results  directly 
determined,  except  for  the  aforementioned  samples  that  may  have  "ink-bottle" 
pores . 

Summarizing,  vacuum  drying  is  a  function  of  the  pore  structure  in 
the  second  falling  rate  period  of  porous  media  and  pore  size  can  be 
calculated  by  equation  (17).  The  advantage  of  this  method  over  other  ways 
of  determining  permeability  is  the  simplicity  of  sample  preparation  and 
the  ease  with  which  the  test  can  be  made.  It  is  easy  to  imagine  automated 
equipment  that  could  be  adapted  from  standard  laboratory  equipment  at  a 
nominal  cost.  The  accuracy  seems  to  be  reasonably  good,  especially  with 
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non-homogeneous  materials  like  aggregates.   In  fact,  one  of  the  signifi- 
cant virtues  of  this  test  method  for  measuring  permeability  is  the  in- 
herent averaging  process  of  the  model;  i.e.,  the  vapor  permeates  radially 
in  two  directions  when  the  specimen  shape  is  cylindrical  as  it  was  in 
this  investigation,  thus  the  resulting  permeability  becomes  an  average  of 
that  in  two  directions.  The  prime  disadvantage  is  the  equation's  depend- 
ence on  tortuosity.   It  is  recommended  that  tortuosity  be  measured  on  the 
same  sample  that  evaporation  tests  are  made  on,  as  was  done  in  this  study. 
This  recommendation  is  extremely  important  owing  to  the  larger  variation 
in  tortuosity  from  specimen  to  specimen,  even  in  those  materials  that 
seem  to  be  relatively  homogeneous. 

Interpretation  and  Application 
It  has  been  recognized  for  some  time  that  the  void  properties  of 
coarse  aggregate  affect  to  a  large  extent  its  performance  in  concrete. 
The  ultimate  scope  of  this  work  was  to  relate  more  accurately  the  function 
of  the  void  system  of  aggregate  with  concrete  performance.  As  demonstrated 
in  the  review  of  the  literature,  no  individual  parameter  defines  categori- 
cally the  potential  performance  of  an  aggregate. 

As  far  as  the  physical  properties  of  pore  size  distribution,  porosity, 
tortuosity,  permeability  and  coefficient  of  diffusion  are  concerned,  a 
significant  correlation  with  the  field  performance  history  and  the  dura- 
bility results  was  lacking.  For  example,  the  rock  designated  as  IS  would 
be  classified  as  a  good  aggregate  if  our  standards  conformed  to  the  cri- 
teria of  low  porosity  and  absorption  but  according  to  field  records  the 
performance  of  IS  has  been  appalling.  These  observations  are  verified 
by  the  results  of  the  ASTM  C291  tests  reported  in  this  investigation. 
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Likewise,  there  are  aggregates  (IDU)  that  have  an  excellent  field  per- 
formance history  and  have  demonstrated  the  same  excellence  in  the  C291 
test.   IDU  has  a  moderately  high  porosity,  a  high  permeability,  relatively 
large  pores,  and  a  pore  size  distribution  that  lies  between  distributions 
of  aggregates  that  have  had  a  very  poor  performance  history.  One  of  the 
significant  conclusions  of  this  study  is  that  the  individual  pore  para- 
meters of  themselves  are  unable  accurately  to  define  the  performance  of 
an  aggregate. 

The  basic  premise  of  freeze-thaw  theory  is  that  the  void  system  must 
be  sufficiently  saturated,  at  least  above  the  critical  point  of  about  91 
percent.  Therefore,  if  an  aggregate  cannot  become  critically  saturated 
it  should  never  cause  freezing  and  thawing  trouble  in  concrete.  Certain 
aggregates  are  unable,  even  under  the  most  severe  laboratory  conditions 
of  vacuum  saturation  to  become  critically  saturated.   Such  results,  shown 
in  Table  6  ,  were  found  in  this  study.  However,  the  fact  still  remains 
that  some  such  aggregates,  such  as  IS,  have  a  poor  performance  history. 

It  is  difficult  to  imagine  a  pore  system  that  would  not  become  filled 
under  vacuum  saturation.  The  total  void  volume  in  this  study  was  measured 
with  the  McLeod  porosimeter,  which  in  essence  pumps  the  air  from  the  void 
system  and  entraps  it  in  a  microburet  where  it  is  measured.  The  question 
is:  why  won't  water  penetrate  into  those  voids  from  which  air  is  readily 
extracted?  This  question  seems  to  merit  further  consideration. 

Owing  partially  to  the  low  permeability  of  the  cement  paste  in  a 
concrete  system  it  has  been  reported  (  4 )  that  the  saturation  time  of  an 
aggregate  in  concrete  may  be  several  months.  To  obtain  an  appreciation 
for  at  least  an  order  of  magnitude  that  may  be  reasonable,  the  following 
hypothetical  problem  was  created  --  How  long  will  it  take  fully  to 
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saturate  a  cube  of  dry  aggregate,  of  1  cnr  volume,  which  has  a  cover  of 
saturated  cement  paste  1  cm  thick  that  is  fully  hydra  ted  and  has  a  water- 
cement  ratio  *  0.5.  The  permeability  of  the  paste  is  assumed  to  be 

-16   2 
1.37  x  10    cm  (59).  This  permeability  is  probably  lower  than  that 

found  in  the  paste  in  some  concrete,  but  it  may  be  realistic  for  pavements. 

Applying  the  above  conditions  to  Darcy's  law,  and  assuming  that  the  driving 

potential  is  developed  by  the  capillarity  of  the  aggregate,  we  have 

ts  =  7.03  x  103e  d  (25) 

where 

ts  «=  time  of  saturation,  hours 
e  =  porosity 

d  =  mean  pore  diameter,  microns 
The  saturation  times,  calculated  from  Eq.  (25)  and  using  the  pore  diameters 
measured  by  evaporation  are  shown  in  Table  13. 

Also  of  interest  is  the  theoretical  drying  time  of  the  aggregate  in 
a  concrete.  The  drying  process  is  considerably  different  than  that  of 
wetting  or  absorption.  First,  sufficient  negative  pressure  must  be 
developed  to  nucleate  (see  Appendix  A)  the  fluid  in  larger  pores.   This 
is  usually  not  difficult.  The  negative  pressure  becomes  primarily  a 
function  of  the  relative  humidity  at  the  air-meniscus  interface,  and  can 

be  evaluated  by  the  Kelvin  equation.  For  example,  a  relative  humidity  of 

6      ,  2 
90  percent  will  develop  139  x  10  dynes/cm  pressure.  This  pressure  must 

be  corrected  by  the  capillary  potential  developed  by  the  aggregate  having 

an  average  diameter  of  d.  Applying  this  criterion  to  the  Darcy  equation, 

we  can  arrive  at  a  theoretical  desaturation  time,  t..,  for  aggregate,  where 

the  subscript  denotes  the  relative  humidity. 

1.45  x  102  ds  (26) 

*90     d  -  0.02  V   ' 
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Table 

13 

Saturation  and 

Drying  Times 

Saturation 

Drying 

Ident. 
No. 

time, 
tg,  hours 

time 
t^j,  hours 

4 

1-1S 

150 

37 

67-2S 

100 

6 

17 

47-2S 

715 

16 

45 

9-1S 

if  200 

30 

140 

MM 

4800 

16 

300 

ML 

2500 

17 

147 

MTS 

1050 

7 

150 

IS 

75 

9 

8 

IA 

146 

8 

18 

II 

1380 

11 

125 

IDU 

3600 

11 

33 

IDE 

150 

5 

30 
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The  drying  times  are  listed  for  each  aggregate  at  90  percent  relative 
humidity  in  Table  13  •  These  values,  though  not  perhaps  exact,  provide  an 
order  of  correctness.  Qualitatively,  these  figures  are  supported  by 
unreported  work,  by  Monfore  (60)  and  Pfeifer  (6l)» 

Shown  along  with  saturation  times  and  drying  times  in  Table  13  are 
the  ratios  of  these.  The  higher  ratio  indicates  that  an  aggregate  remains 
in  an  uncritical  state  of  saturation  for  longer  periods.  However,  it 
appears  that  for  aggregates  having  ratios  of  up  to  300,  their  saturation 
may  be  enough  to  cause  frost  damage.  One  can  envisage  an  aggregate  with  a 
void  system  whose  pore  area  is  large  enough  to  be  immune  from  the  point  of 
critical  saturation.  This  state  would  be  reflected  in  the  saturation- 
drying  ratio.  Some  lightweight  aggregate  may  be  in  this  category.  Also, 
it  is  believed  that  perhaps  this  type  of  protection  is  offered  to  1-1S 
and  IDU,  even  though  the  data  and  results  do  not  clearly  indicate  it. 
This  may  be  explained  partially  by  the  inability  of  the  mercury  porosimetry 
equipment  to  measure  the  macro-pores.  They  definitely  exist,  since  they 
can  be  seen  by  the  unaided  eye. 

The  above  discussion  has  indicated  that  most  carbonate  aggregates 
have  void  systems  that  are  capable  of  acquiring  water  at  a  rate  sufficient 
to  remain  oritically  saturated  under  several  conditions.  Aggregates  that 
have  been  classified  as  good  and  bad  are  equally  capable  of  becoming 
saturated  under  similar  circumstances. 

Thus  it  behooves  us  to  demonstrate  why  some  aggregates  are  able  to 
protect  themselves  against  freezing  and  thawing,  while  others  are  not. 
When  fully  saturated  aggregates  are  frozen,  the  aggregate-paste  system 
must  accommodate  the  increase  of  volume,  V  , 

V8  ■  0.09 e  (27) 
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per  unit  bulk  volume, or  the  system  will  destroy  itself.  Three  processes 
have  been  postulated.  These  mechanisms  are  usually  referred  to  as: 
(1)  elastic  accommodation,  (2)  critical  size,  and  (3)  expulsion  into  paste. 
Much  of  the  following  discussion  is  after  Verbeck  and  Landgren  (4). 

In  general  those  aggregates  that  have  low  porosities  and  permea- 
bilities have  shown  good  performance.  An  aggregate  is  able  to  accommodate 
elastically  some  of  the  excess  volume  when  the  water  in  its  void  system 
is  frozen.  The  elastic  accommodating  volume,  V  ,  depends  to  a  large 
extent  on  the  tensile  strength,  P. ,  of  the  aggregate  and  on  its  modulus 
of  elasticity,  E.  The  accommodating  volume  per  unit  bulk  volume  can  be 
estimated  by 

3  Pt  (1  -  2/0 
VA  =  — ±~ g (28) 

These  "volumes",  which,  strictly,  are  not  volumes  but,  like  V  ,  are  ratios 
of  volumes  and  dimensionless,  are  listed  in  Table  14  for  each  aggregate. 
The  tensile  strengths,  F  ,  measured  by  the  tensile  splitting  test,  are 
given  in  Table  6  •  The  moduli  of  elasticity  and  Poisson*s  ratios  were 
assumed  to  be  3  x  10  psi  and  1/4  respectively.  If  V  <  V.  the  aggregate 

©  A 

is  able  elastically  to  accommodate  the  excess  volume.  This  comparison, 
shown  in  Table  14  for  the  aggregates  investigated,  indicates  that  none  of 
these  aggregates  is  able  elastically  to  protect  itself.  Therefore,  we 
must  assume  that  some  other  process  has  aided  in  the  protection  of  those 
aggregates  that  are  classified  as  good. 

If  a  saturated  aggregate  is  unable  internally  to  accommodate  the 
volume  change  developed  during  freezing,  two  outcomes  are  possible, 
depending  upon  the  freezing  rate  and  the  aggregate's  permeability. 
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By  using  Darcy's  law  a  "critical  size"  of  an  aggregate  can  be  esti- 
mated as  follows  (1*): 

27.7  K  P, 


(29) 


max   0.09  €  dF/dt 

where : 

L     =  maximum  permissible  size,  in. 
max 

dF/dt  =  freezing  rate,  cm/sec 

P     =  tensile  strength  of  aggregate,  psi 

K     =  permeability,  cm/ sec 
The  critical  size  calculated  for  each  aggregate  is  shown  in  Table  14.  The 
permeabilities  used  in  the  calculation  are  those  from  Table  10  .  The 
freezing  rate,  taken  from  Verbeck  and  Landgren  (4),  was  1.5  inches  per 
1  hr.  (1.06  x  10   cm/sec).  They  considered  this  rate  to  approximate 
those  encountered  in  natural  exposures  and  in  some  laboratory  tests. 
Gress  (62)  measured  a  freezing  rate  of  2.2  in./hr.  for  3  x  U  x  16  inch 
prisms  tested  in  the  equipment  used  for  the  freeze-thaw  tests  reported 
in  this  study. 

This  critical  size  is  the  maximum  distance  water  can  be  forced 
through  such  a  material  by  the  advancing  ice  front.  It  is,  therefore, 
the  maximum  size  the  aggregate  piece  can  be  without  fracturing  internally 
when  frozen.  Aggregates  that  typically  have  such  a  short  critical 
distance  that  it  falls  within  usual  particle  sizes  of  coarse  aggregate 
are  the  lightweight  cherts  and  shales  that  so  frequently  exhibit  pop-out 
behavior. 

From  the  results  in  Table  14 ,  it  is  evident  that  aggregate  IS  has 

a  much  shorter  critical  distance  than  any  other.  This  finding  is  signif- 
icant. From  the  freeze-thaw  results  in  Fig.  E5it  is  observed  that 
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particles  of  IS  have  fractured  internally  into  relatively  small  pieces. 
This  failure  was  unique  among  the  specimens  studied. 

An  obvious  question  is,  why  does  this  aggregate  fail  to  have  pop 
outs,  which  are  usually  associated  with  small  critical  size  and  internal 
failure?  It  is  the  opinion  of  this  writer  that  the  surrounding  paste  is 
able  to  accommodate  elastically  the  increased  volume  change  enough  to 
protect  the  paste-aggregate  unit  from  forming  the  notable  cone  fracture 
that  is  usually  associated  with  this  failure.  Most  pop-out  failures  due 
to  freezing  and  thawing  are  associated  with  aggregates  that  have  low 
permeabilities  and  high  porosities,  and  that  are  able  to  generate  enough 
force  that  the  paste  Is  as  unable  as  the  aggregate  to  accommodate  the 
volume  change  and  if  they  are  near  the  surface,  a  cone  fracture  develops. 

The  other  aggregates  are  able  to  expel  their  excessive  water  on 
freezing  into  the  surrounding  paste  without  internal  fraoture.  Therefore, 
the  paste  must  be  capable  of  accommodating  this  expelled  water  or  the 
hydraulic  pressure  developed  in  the  paste  will  rupture  it  close  to  the 
interface  in  what  is  essentially  a  bond  failure. 

The  "expulsion  into  the  paste"  phenomenon  becomes  primarily  a 
question  of  freezing  rate,  aggregate  size,  and  porosity,  and  the  permea- 
bility and  air  content  of  the  surrounding  paste.  Such  a  process  can  best 
be  analyzed  by  assuming  a  spherical  saturated  aggregate  covered  with  an 
air-entrained  paste. 

The  spherical  shell  of  paste  can  accommodate  the  following  volume  of 
expelled  water  from  the  aggregate 

V  o  A  4  tt  r2t  -  0.1  e  4  n  A        (30) 

c 
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where 

V  =  volume  of  space  in  paste  available  "to  expelled  water 

A  «=  fractional  air  content  of  the  paste 

r  =  radius  of  the  spherical  aggregate  particle 

t  m     thickness  of  the  paste  shell 

e  =  capillary  porosity  of  the  paste 

If  this  value  is  then  divided  by  the  volume  of  the  aggregate  sphere 
a  dimensionless  parameter,  V™,  is  obtained 

*-# 

that  can  be  compared  with  the  values  of  V  of  water  expelled  from  unit 
volume  of  the  aggregate. 

VT  was  determined  by  using  the  following  values:  A  =  0.18,  e  =  0.15, 
r  =  1  cm,  and  t  =  0.023  «n»  which  is  the  value  that  has  been  assumed 
reasonable  for  the  critical  distance  for  water  flow  in  a  paste  during 
freezing.  The  value  obtained  is  11.3  x  10  .   This  value  is  compared  with 
the  individual  values  of  V.  in  Table  14  .  If  V-  is  larger  than  V  an 
aggregate  piece  of  this  size  (r  =  1  cm)  should  be  protected  from  damage 
due  to  expulsion  of  water  into  the  surrounding  paste,  and  vice  versa. 

For  the  preceding  discussion  the  mechanisms  are  seen  to  rely  on 
several  assumptions.  The  inferences  generated  from  these  concepts  are 
probably  only  semi-quantitative.  The  results,  which  are  summarized  in 
Table  14,  do  delineate  accurately  the  behavior  of  all  the  aggregates 
with  the  sole  exception  of  MTS.  Some  reports  have  shown  that  the  perform- 
ance of  MTS  has  been  acceptable  in  some  pavements,  but  it  is  generally 
considered  a  poor  material.  The  accelerated  laboratory  tests  show  an 
expansion  factor  of  7.0  and  a  durability  factor  of  40.  Acoording  to  these 
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values,  this  aggregate  would  be  classified  marginal  to  poor.  It  is 
speculated  that  there  is  wide  variance  in  the  properties  of  MTS,  sufficient 
that  some  of  its  geological  formations  may  be  acceptable  while  other 
portions  may  not. 

Hie  results  of  the  accelerated  freezing  and  thawing  tests,  ASTM 
Method:  C291,  are  of  interest  when  compared  with  the  field  service 
records  of  the  materials.  When  each  aggregate  was  received  from  its 
source,  a  life  expectancy  for  highway  pavements  was  ascertained  a  priori. 
The  life  expectancy,  for  all  aggregates  exoept  47-2S  and  67-2S  is  plotted 
against  the  expansion  factor  in  Fig.  30.  This  correlation  is  good.  Some 
significance  can  perhaps  thus  be  attached  to  the  expansion  factor  as  a 
quantitative  parameter  for  expressing  durability. 

Aggregates  in  this  study  that  have  expansion  factors  greater  than  8j0 
have  had  good  field  performance.  Values  between  7.0  and  8j0  are  marginal. 
Aggregates  with  expansion  factors  lower  than  7.0  have  shown  poor  perform- 
ance. The  durability  factor,  which  is  the  most  commonly  used  index  of 
frost  susceptibility,  did  not  delineate  the  performance  of  the  aggregate 
as  well  as  the  expansion  factor  did.  This  observation  is  apparent  when 
the  expansion  data  and  the  dynamic  modulus  of  elasticity  data  are  compared 
(Figs.  26  to  29  and  Figs.  E6  to  E8). 

The  photographs  of  the  freeze-thaw  specimens  portray  some  of  the 
apparent  failure  mechanisms  of  these  aggregates  in  concrete.  Listed  with 
the  results  of  the  freezing  and  thawing  data  is  a  qualitative  assessment 
of  the  failure  type.  These  failures  were  classified  as  either  boundary 
distress  or  aggregate  internal  fracture  (Table  12). 

The  first  assessment  of  the  boundary  distress  would  classify  it  as 
a  bond  failure  near  the  aggregate  and  paste  interface,  but  on  close 
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examination  of  the  failure  surface,  particularly  on  the  specimen  designa- 
ted as  9-lS(C),  a  thin  shell  of  aggregate  had  adhered  to  the  paste  surface. 
This  type  of  phenomenon  has  been  termed  in  Iowa  as  an  "onion  skin  failure. * 
The  surface  of  failure  seems  to  lie  a  small  distance  into  the  aggregate* 
It  has  been  suggested  (63 )  that  while  the  concrete  is  still  plastic ,  the 
paste  is  able  to  migrate  a  short  but  significant  distance  into  the  aggre- 
gate void  system,  and  that  the  point  of  penetration  becomes  an  incipient 
failure  surface.  Undoubtedly  the  permeability  of  this  paste  zone  in  the 
aggregate  is  reduced.  If  the  pore  diameter  of  the  aggregate  void  system 
lies  in  the  1  to  20  micron  range,  a  filtering  process  can  be  envisaged. 
A  process  where  the  small  cement  particles  and  soluble  hydration  products 
are  permitted  to  filter  into  the  aggregates'  void  system,  but  the  entrained 
air  bubbles  are  filtered  out.  It  is  possible  that  this  zone  will  have 
permeability  equivalent  to  the  paste  system  but  be  unprotected  by  air- 
entrainment,  and  thus  allowing  the  zone  to  become  more  vulnerable  to  the 
expelled  water  which  cannot  be  accommodated  in  the  zone. 

The  internal  fracture  failure  was  evident  only  in  aggregate  IS.  The 
remaining  failures  appeared  to  be  the  boundary  distress  types.  These 
observations  are  in  accordance  to  the  predicted  mechanisms  listed  in 
Table  14  . 
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SUMMARY  OF  RESULTS 

The  major  findings  of  this  study  can  be  summarized  as  follows: 

1.  The  technique  used  to  determine  the  rate  of  evaporation  of  a 
liquid  from  the  pores  of  a  porous  solid  was  simple,  and  the  data  were 
reproducible. 

2.  An  empirical  equation  was  written  that  describes  the  results  of 
the  air-drying  experiments  up  to  a  relatively  dry  state  of  the  sample. 

A  formal  diffusion  coefficient  was  calculated  from  these  results,  but  no 
relationship  with  the  pore  structure  of  the  solid  was  found. 

3.  The  air-drying  experiments  with  the  thin,  fritted-glass  samples 
showed  a  long,  constant-rate  period,  the  only  found  for  any  samples.  The 
evaporation  rate  was  approximately  the  same  for  samples  the  pore  size  of 
which  varied  by  an  order  of  magnitude. 

k.     When  the  results  of  the  vacuum  evaporation  experiments  were 
interpreted  by  means  of  a  modified  slip-flow  equation  the  values  of  pore 
size  obtained  were  frequently  about  the  same  as  those  obtained  by  the 
mercury  intrusion  experiments  or  by  a  direct  measurement  of  permeability. 

5.  Ho  correlation  was  found  between  any  single  pore  parameter  of 
the  aggregates  and  their  performances  in  laboratory  freezing  and  thawing 
tests  of  concrete  or  their  field  histories. 

6.  The  field  performance  of  the  aggregates  was  reasonably  assessed 
by  laboratory  weathering  tests.   Poor  aggregates  were  characterized  by  a 
loss  in  dynamic  modulus  of  elasticity,  a  length  increase,  and  a  weight 
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increase  of  the  concrete  samples  containing  them.  The  weight  gain  was 
probably  due  to  an  increased  absorption  of  water  resulting  from  the 
development  of  cracks  in  the  concrete.   It  was  indicated  that  length 
change  was  an  equally  good,  if  not  better  parameter  than  change  in 
dynamic  modulus  (and  the  resulting  durability  factor)  as  an  interpretation 
of  the  laboratory  freezing  and  thawing  tests. 

7.  By  an  application  of  the  data  developed  in  this  study  and  reason- 
able assumptions,  the  field  performances  of  all  but  one  of  the  aggregates 
could  be  characterized  in  terms  of  the  mechanisms  set  forth  by  Verbeck 
and  Land gr en  (k) . 

8.  It  was  indicated  that  aggregates  in  concrete  can  become  highly 
saturated  fairly  readily. 
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CONCLUSIONS 

Based  on  the  results  of  this  study,  the  following  conclusions  seem 
reasonable.  They  cannot  be  expected  necessarily  to  apply  to  types  of 
materials  other  than  those  tested. 

1.  The  steady-state,  slip-flow  equation,  (G7) ,    can  be  modified  by 
the  inclusion  of  another  porosity  term  to  give  equation  (17) >   that  is 
applicable  to  describe  the  vacuum  evaporation  process  in  the  second  drying 
rate  period. 

2.  Vacuum  evaporation  of  water  from  a  porous  solid  is  a  simple 
technique  that  can  be  used  to  determine  values  of  pore  size  and  permea- 
bility that  are  comparable  to  those  determined  by  mercury  intrusion  or 
direct  measurement. 

3.  The  constant  rate  of  evaporation  during  air  drying  is  not 
related  to  the  pore  size  of  the  porous  material. 

k.  No  single  pore  parameter,  including  the  pore  size  distribution 
function,  of  a  carbonate  aggregate  is  a  good  indication  of  its  behavior 
in  freezing  and  thawing  exposure  in  concrete.  , 

5.  The  mechanisms  of  Verbeck  and  Landgren  (k)   can  be  used  with 
easily  obtained  laboratory  test  results  accurately  to  characterize  the 
probable  field  performance  of  carbonate  aggregates  in  freezing  exposure. 
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APPENDIX  A 

Influence  of  Ion  Adsorption  on  Tortuosity  Measurements 
for  Materials  Having  High  Internal  Surface  Areas 

Materials  whose  internal  surface  areas  are  large  usually  yield  high 
tortuosity  values  when  evaluated  by  electrical  resistivity  methods.  For 
example,  Vycor  has  a  measured  tortuosity,  K^,  that  exceeds  1000.  Such 
values  are  not  realistic.  If  one  interprets  this  value  literally,  the 
tortuous  path,  or  the  length  of  fluid  travel,  would  be  about  thirty 
times  the  length  of  the  sample. 

When  an  electrical  analogue  is  used  for  tortuosity  measurements,  the 
resistivity  of  the  conducting  fluid,  5  percent  sodium  chloride  solution, 
is  assumed  to  be  the  same  in  the  pore  system  as  in  the  bulk  solution, 
and  the  solid  material  is  assumed  to  be  non-conducting.   The  following 
calculation  shows  that  this  assumption  may  be  incorrect,  especially  for 
materials  having  large  internal  surface  areas. 

Assume  the  sample  to  consist  of  one  gram  of  Vycor  porous  glass. 
The  porosity  will  be  about  0.35  and  the  specific  surface  will  be  about 
250  m  /g.  Both  values  are  from  the  manufacturer.  When  the  sample  is 
saturated  with  the  conducting  solution,  5  percent  NaCl,  assume  the  ions 
are  adsorbed  on  the  solid  surface  and  that  equal  amounts  of  cation  and 

anion  are  sorbed.   If  the  ions  are  simply  sorbed,  i.e.  without  any 

2 
hydration  hull,  they  will  occupy  perhaps  10  A  ,  on  the  average.  A 

21 
complete  monolayer  of  ions  will  then  involve  the  abstraction  of  2.5  x  10 
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ions.  If  one  assumes  a  bulk  density  of  1*5   g/cnr  for  the  glass,  the  one- 

3 
gram  sample  will  have  a  pore  volume  of  0.233  cm  .  If  this  pore  volume 

is  saturated  with  a  5  percent  NaCl  solution  the  pores  will  contain 

20 
initially  about  2.5  x  10   ions.   In  other  words,  if  the  above  is  a 

reasonable  proposition,  the  surface  of  such  a  high-area  solid  has  the 

ability  to  abstract  more  than  all  the  ions  added.  Such  an  event  would 

result  in  a  great  lowering  of  the  conductivity,  which  would  be  interpreted 

by  the  conventional  calculation  as  a  very  high  tortuosity. 

Ihe  question  of  the  tortuosity  of  Vycor  was  investigated  briefly 

by  means  of  absorptivity  measurements.  Absorptivity  was  defined  by 

Powers  and  Brownyard  ( 6*0  as 

f  -   K.t  (Al) 

where 

V  =  volume  of  a  liquid  absorbed  by  a  regularly-shaped 
specimen 

A  =  cross-sectional  area  normal  to  the  liquid's  flow 

t  =  time  of  absorption 

Kfi  =  absorptivity  coefficient 

Such  a  measurement  was  made  on  the  Vycor  samples  using  the  method 
of  Dolch  (21).   Ihe  results  are  given  in  Fig.  Al  .  From  the  data  in  the 
figure,  V/A  =  0.024  at  t  =  1  min. 

K  =  iV^i  =  {°'lf)2   =  9.6  x  10"6  cm2/sec      (A2) 

Dolch  ( 2l)  showed  the  following  relationship  between  absorptivity 
and  tortuosity 
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where    K   =  Tortuosity 

Y      =  Surface  tension 
e      =  Porosity 
I     =  Viscosity 
Solution  of  the  equation  gives  a  K.  of  10.6.   This  would  yield  a 
tortuous  path  length  of  3*2  or  10  times  less  than  the  value  measured  by 
the  electrical  analogue. 
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APPENDIX  B 

Ihe  Mathematical  Relationship  Between  the  Mass-rate 
Function  and  Time-rate  Function 

Drying  parameters  are  commonly  expressed  as  the  "Rate"  of  evapora- 
tion versus  "Percent  Moisture*1  content  or  "The  Mass-rate  Function. "  As 
shown  in  Fig.  9  of  the  text,  the  first  falling  rate  period  has  a  para- 
bolic curve.  If  "Time"  and  "Rate"  of  evaporation  are  plotted  on  a  log- 
log  scale,  the  plot  of  the  first  falling  rate  period  becomes  a  straight 
line,  and  from  this  simple  relationship  the  following  equation  can  be 
written: 

V=Ktb=Va  C*1) 

A         a 

where 

t  =  time 

V  =  volume  of  fluid  lost  by  evaporation 

A  =  exposed  area  of  sample 
The  rate  of  evaporation,  R,  per  unit  area  can  be  found  by  differentiating 
equation  1  with  respect  to  t 

f!a  =  RBKbtW  (B2) 

dt 

By  eliminating  t  from  equations  1  and  2 

R  =  CV^r1  (B3) 

D 

All  values  in  equation  B3can  be  readily  obtained  from  the  Time-rate 
function  thus  permitting  the  parabolic  curve  generated  from  the  Mass- 
rate  function  to  be  expressed  mathematically. 


APPENDIX  C 


NUCLEATION  RADII 


103 


APPENDIX  C 


Nucleation  Radii 


Most  of  the  porous  systems  in  concrete  are  not  analogous  to  the 
"bundle  of  capillaries"  concept.  A  more  likely  model  would  be  a  system 
of  isolated  cavities  inter-connected  by  a  series  of  smaller  capillaries. 
Going  on  the  assumption  that  fluid  evaporates  only  from  the  surface  of  the 
capillary  meniscus,  none  of  the  water  from  the  isolated  cavities  would  be 
available  for  removal  until  the  evaporating  front  had  retreated  to  the 
vicinity  of  the  cavity. 

According  to  Powers  (65),  several  data  show  that  these  larger  isolated 
cavities  pre-empt  the  smaller  capillaries.  Powers  used  a  theory  developed 
by  Bernath  (66)  to  show  a  possible  model  by  which  these  larger  isolated 
cavities  could  be  emptied  while  they  were  yet  surrounded  by  capillaries 
completely  full  of  water. 

From  the  perfect  gas  equation  of  state,  for  a  bubble  of  vapor  in 
such  a  cavity 

Pb  (|  sr3)  -  n^kT  (C-l) 

where  T  is  the  absolute  temperature ;  P.  represents  the  gas  pressure  in 
the  cavity  that  is  assumed  to  be  a  sphere;  k  is  Boltzman's  constant  and 
n.  is  the  number  of  molecules  in  the  bubble.  The  cavity  pressure  is  also 
equivalent  to  the  external  pressure, 

PK  "  *r   +  ¥  (C-2) 

b    c   r 


1C& 


where  P  is  the  capillary  hydrostatic  pressure,  in  tension,  and  Y    is  the 


surface  tension  of  the  fluid.  Substituting  from  equation  C-l  we  have 

(C-3) 


(f+  V  C^nr^-VcT 


C  ^3 

This  equation  represents  the  free  energy  content  of  the  cavity  at  a 
given  temperature  as  a  function  of  the  radius,  r.  This  function  is  such 
that  there  exists  a  critical  value  of  the  radius,  r',  which  can  be  obtained 
by  differentiating  equation  C-3  and  equating  it  to  aero. 


...*  -* 


3  Pe 


(CJO 


The  equation  of  equilibrium  is : 


2  Y 


ZY 


rb     c  c 


<< 


(C-5) 


Let  us  illustrate  this  relationship  by  an  example.  If  drying  produces 
a  relative  humidity  of  92  percent,  a  meniscus  having  a  radius  of  curvature 
of  12^A  will  be  developed  (Kelvin  equation)  which  in  turn  will  generate 
in  the  pore  system  a  negative  pressure  of  116  atm  which  is  sufficient  to 
nucleate  radii  larger  than  89A  ,  allowing  5A  for  adsorbed  film. 
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Table  Dl  -  Drying  Data  for  Vacuum  Evaporation 

Specimen  Identification 

67-2S         U7-2S         9-1S  MM(1)         MM(2) 


AT,* 

AW,* 

AT, 

AW, 

AT, 

AW, 

AT, 

AW, 

AT, 

AW, 

min 

gram 

min 

gram 

min 

gram 

min 

gram 

min 

gram 

4.0 

0.080 

k.o 

0.125 

4.0 

0.321 

4.0 

0.219 

4.0 

0.347 

4.o 

0.010 

k.o 

0.043 

4.0 

0.222 

4.0 

0.086 

4.0 

0.190 

It-.O 

0.005 

k.o 

0.030 

4.0 

0.142 

4.0 

0.0U8 

4.0 

0.094 

4.0 

0.005 

k.o 

0.025 

4.0 

0.103 

4.0 

0.030 

4.0 

0.050 

k.o 

0.004 

k.o 

0.020 

4.0 

0.082 

4.0 

0.028 

4.0 

0.035 

k.o 

0.00^ 

k.o 

0.020 

4.0 

0.073 

4.0 

0.017 

10.0 

0.057 

k.o 

0.003 

k.o 

0.015 

4.0 

0.055 

4.0 

0.013 

10.0 

0.029 

k.o 

0.003 

6.0 

0.021 

4.0 

0.046 

10.0 

0.021 

10.0 

0.015 

8.0 

0.003 

6.0 

0.019 

8.0 

0.077 

32.0 

0.027 

20.0 

0.004 

22.0 

0.009 

8.0 

0.019 

10.0 

0.060 

15.0 

0.006 

70.0 

0.005 

32.0 

0.007 

12.0 

0.023 

10.0 

0.052 

15.0 

0.003 

46.0 

0.007 

12.0 

0.018 

20.0 

0.047 

k6.o 

O.OO6 

12.0 

0.016 

60.0 

0.021 

36.0 

0.004 

12.0 

0.0l3 

4o.o 

0.001 

12.0 

12.0 
20.0 
20.0 
20.0 

0.011 

0.009 
0.013 
0.011 
0.009 

*  AT  is  the  time  interval  in  minutes  in  which  a  weight  change  of  AW  in 
gram  was  lost. 


106 


Table  Dl  (Cont'd.) 
Specimen  Identification 
ML(1)        ML(2)  ML(3)        M.(4)        MTS(l) 


AT, 

AW, 

AT, 

AW, 

AT, 

AW, 

AT, 

AW, 

At, 

AW, 

min 

gram 

min 

gram 

min 

gram 

min 

gram 

min 

gram 

4.0 
4.0 
4.0 
k.O 

h.o 

0.238 
0.100 
0.071 
0.063 
0.040 

k.O 
k.O 
k.O 
k.O 
k.O 

0.225 
O.O85 
0.052 
0.038 
0.029 

4.0 
4.0 
4.0 
4.0 
4.0 

0.282 
0.107 
0.059 
0.040 

0.037 

4.0 
4.0 
4.0 
4.0 
4.0 

0.232 
0.192 
0.082 
0.042 
0.027 

4,0 
4.0 
4.0 
4.0 
4.0 

0.108 
0.034 
0.024 
0.014 
0.014 

5-0 
3.0 
k.O 
8.0 
8.0 

o.o4i 

0.023 
0.024 
0.036 
0.026 

k.O 
k.O 
k.O 
k.O 
k.O 

0.025 
0.020 
0.018 
0.014 
0.014 

4.0 
4.0 
4.0 
5.0 
13.0 

0.023 

0.021 

0.019 
0.015 
0.027 

4.0 
4.0 
5.0 
7.0 
43.O 

0.027 
0.027 
0.023 
0.021 
0.056 

4.0 
4.0 
4.0 
8.0 
8.0 

0.010 
0.007 
0.005 
0.010 

0.006 

12.0 
12.0 
12.0 
12.0 
12.0 

0.030 
0.023 
0.013 
0.011 
0.010 

8.0 
12.0 
12.0 
12.0 
12.0 

0.0  is 
0.023 
0.015 
o.oi4 
o.oli 

25.0 
25.O 
50.0 

0.032 
0.019 

0.020 

0.018 
0.005 

12.0 
12.0 
12.0 
24.0 
24,0 

0.007 
0.007 
0.004 

0.007 
0.004 

12.0 
12.0 

0.006 
0.005 

12*0 
12.0 
12.0 

0.009 
0.006 
0.006 

12.0 

0.003 
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Table  Dl  (Cont'd.) 
Specimen  Identification 
MTS(2)  IS  11(1)  11(2)         IDU(l) 


AT, 

AW, 

AT, 

AW, 

AT, 

AW, 

AT, 

AW, 

AT, 

AW, 

min 

Sra"» 

min 

gram 

min 

gram 

min 

gram 

min 

gram 

4.o 

0.106 

4.0 

0.024 

4.0 

0.196 

4.0 

0.125 

4.0 

0.132 

4.0 

0.032 

k.o 

0.006 

4.0 

O.O69 

4.0 

0.038 

4.0 

0.047 

k.o 

0.023 

k.o 

0.005 

4.0 

0.043 

4.0 

0.026 

4.0 

0.026 

4.0 

0.020 

k.o 

0.004 

4.0 

0.031 

4.0 

0.020 

4.0 

0.019 

k.o 

0.015 

8.0 

0.007 

4.0 

0.025 

4.0 

0.013 

4.0 

0.014 

k.o 

0.012 

12.0 

0.005 

4.0 

0.019 

4.0 

0.010 

4.0 

0.012 

k.o 

0.008 

12.0 

0.005 

4.0 

0.015 

4.0 

0.010 

4.0 

0.008 

k.o 

0.007 

22.0 

0.007 

4.0 

0.011 

4.0 

0.006 

4.0 

0.007 

8.0 

0.011 

50.0 

0.013 

8.0 

0.022 

8.0 

0.010 

4.0 

0.006 

8.0 

0.007 

62.0 

0.014 

8.0 

0.010 

8.0 

0.008 

4.0 

0.005 

12.0 

0.007 

124.0 

0.018 

14.0 

0.017 

14.0 

0.010 

8.0 

0.007 

12.0 

0.007 

720.0 

0.054 

14.0 

0.009 

14.0 

0.006 

14.0 

0.010 

14.0 

a.  007 

24.0 

0.015 

24.0 

0.010 

14.0 

0.007 

10.0 

0.003 

30.0 

0.100 

30.0 

0.010 

64.  0 

0.015 

30.0 

0.08 

24.0 

30.0 
30.0 

0.005 

0.009 
0.003 
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Table  Dl   (Cont'd.) 


Specimen 

Identification 

IDU(2) 

IDE 

IA( 

10 

IA 

(2) 

AT, 

AW, 

AT, 

AW, 

AT, 

AW, 

AT, 

AW, 

min 

gram 

mtn 

gram 

min 

gram 

min 

gram 

4.0 
4.0 
4.0 
k.O 
4.0 

0.254 

0.103 
0.067 
0.032 
0.020 

k.O 
k.O 
k.O 
k.O 
k.O 

0.057 
0.019 
0.013 
0.010 
0.008 

4.0 
4.0 
4.0 
4.0 
8.0 

0.045 
0.016 
0.010 
0.009 
0.013 

4.0 

4.0 
4.0 
4.0 
4.0 

0.048 

0.017 
0.012 

0.009 
0.009 

k.O 
k.O 
k.O 

k.O 
k.O 

0.011 
0.009 
0.004 

0.003 

0.004 

k.O 
k.O 
k.O 
8.0 
8.0 

0.007 
0.006 
0.005 
0.007 
0.005 

8.0 
18.0 
17-0 
35-0 
25.O 

0.008 
0.014 
0.009 
0.015 
0.006 

4.0 
4.0 
4.0 
8.0 
8.0 

0.008 
0.007 
0.002 
0.008 
0.006 

8.0 

0.000 

14.0 
12.0 

56.0 
43.0 

0.008 
0.005 
0.018 
0.008 

0.042 

17.0 
35-0 
25.0 
45.0 
45.0 

0.007 
0.012 
0.005 
0.007 
0.004 

Fig.  El  Freezing  and  Thawing  Specimen,  MTS 


Fig.  t.2     Freezing  and   Thawing  Specimen,  KL 


Fig.  E3  Freezing  and  Thawing  Specimen,  9-1S 
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Fig.  E5  Freezing  and  Thawing  Specimen,  IS 
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APPENDIX  F 

Geometry  Considerations 

The  purpose  of  this  section  is  to  consider  what  effect  specimen 
geometry  has  on  evaporating  phenomenon.  Cylindrical  specimens  were  used 
almost  exclusively  throughout  this  study.   The  cylindrical  shape  was 
selected  for  two  reasons:  1)  it  could  be  fabricated  and  sized  with  a 
minimum  of  effort  by  coring  techniques,  and  2)  the  cylindrical  shape 
provided  an  optimum  surface  area  for  evaporation,  while  allowing  reason- 
able mathematical  treatment  of  fluid  flow. 

Equation  (1?)  the  vapor  flow  is  assumed  to  move  unidirectionally  and 
perpendicular  to  the  plane  of  evaporation  area  of  which  remains  constant 
throughout  the  process.  However,  with  vapor  flow  from  a  cylindrical 
surface,  the  flow  direction  is  radial  and  normal  to  the  evaporating 
surface,  the  area  of  which  continues  to  diminish.  It  will  be  shown  that 
these  two  processes  are  related. 

Referring  to  Fig.  Fl 

A  =  2  n  \  r  (F-l) 

where  A  is  the  area  of  the  evaporating  surface  at  r.  If  A  is  considered 
to  be  constant,  as  is  the  case  for  the  value  of  "b"  shown  in  Table  8 

VxE2nU  <F-2 

Referring  once  again  to  Fig.  Fl,  and  defining 

A^  =  tt  (a  -  r  ) 

2 

A2-nr 
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FIG.     Fl        SCHEMATIC    OF   EVAPORATION     FROM    A 
CYLINDRICAL     SPECIMEN 
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we  can  write 


2    2      2 
n  (a  -  r  )  m  tt  r  (F-3) 


by  assuming  that  the  average  effect  of  this  process  can  be  represented  at 
the  midpoint  of  volume  loss.  From  expression  (F-3) 


r  =&  CM) 


or  the  effective  area,  A,  becomes 

-       2  n    X      a  .      „. 

A  =  ~~7z (F-5) 

The  ratio  between  Amax  (F-2)  and  A  (F-5)  is  the  correction  factor  -/2. 
This  factor  was  evaluated  experimentally  and  was  found  to  lie  between 
1.40  and  1.50. 

Thus,  the  values  of  "bn  in  Table  8  were  increased  by  the  V2  to 
correct  for  geometry  effects  in  the  pore  size  calculations. 
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APPENDIX  G 

Flow  Equation 

Carman  (18)  recognizes  three  general  regimes  when  a  gas  is  flowing 
through  a  tube  or  porous  medium.  These  are  l)  a  purely  viscous  flow 
region  to  which  the  Poiseuille  equation  applies  and  the  velocity  of  the 
gas  at  the  tube  boundary  is  zero,  2)  a  "slip  flow"  region  in  which  the 
velocity  at  the  tube  wall  is  not  zero  and  which  occurs  when  the  mean  free 
path  of  the  gas  is  about  the  size  of  the  channel  through  which  it  is 
flowing,  and  3)  the  Knudsen  flow  region  in  which  all  viscous  effects 
disappear  and  which  occurs  when  the  mean  free  path  of  the  gas  is  large 
compared  with  the  size  of  the  tube. 

Assume  a  porous  medium  of  approach  area  A  normal  to  the  incoming 
fluid  flow  path  and  length  L  between  inlet  and  outlet  faces.  The  pore 
system  consists  of  n  pores  of  non-circular  cross  section,  the  hydraulic 
radius  of  which  is  m.  The  pores  have  a  real,  or  tortuous,  length  L  which 
is  the  overall  distance  the  fluid  traverses  through  the  medium. 

Carman's  (l8)  equation  for  flow  through  one  such  pore  is 

*1P1   m2?    ,  *  glf  ,„* 

*     "Vle   ^e 

where 

li-  ■  flow  velocity  measured  at  pressure  P. 

AP  ■  pressure  differential  between  inlet  and  outlet  faces 

P  «*  mean  pressure  across  sample 
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K  ■  shape  factor  for  viscous  flow 
,  o 

V    ■  viscosity  of  fluid 

3*  2"fl 
8  »  *V  (—J — )  where  f   is  the  fraction  of  molecules  diffusely 

reflected  when  they  hit  the  pore  wall 

-  /8rt 

v  =  mean  thermal  speed  of  molecules,  i.e.    -rr~ 

KM 


and 


K.  ■  shape  factor  for  slip  component 

Since  the  parameter  measured  in  this  work  is  weight  change, the  volume 
flow  rate  must  be  converted  to  a  rate  of  weight  change  by 

dV 

"i  ■  i  •  -i  (o2> 

where   a  "   cross  sectional  area  of  tube 

V.  ■  volume  of  gas  flowing  through  tube,  measured  at  P1 
and 

t  =  time 

Then,  if  the  gas  is  assumed  to  be  perfect 

W  RT 
plvl  "  ~  <Q3) 

where 

W  -  weight 

M  ■'  molecular  weight  of  gas 

R  =  gas  constant 
and 

T  ■  absolute  temperature 

If  (G3)  is  differentiated  and  combined  with  (G2)  the  result  is 

dV1    ^1   RT 

-d7--dc-<f^-^  («*> 
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The  total  mass  flow  rate  through  the  n  pores  will  be 

dW   ndw   n^laPlM 


dt     dt     RT 


(05) 


or 


2-        k   8  m 
1_  dW  m     naM  ,mfp__        1   ,  („ 

AP  dt      RT  kKqV   Le     3K.Le  ;  K^° ' 

If  the  substitution  of  Ac  ■»  na,  where  e  is  the  porosity,  is  made  and 

o 
the  definition  of  tortuosity,  K  -  (L  /L)  ,  is  inserted,  equation  (G6) 

becomes 

*  -  JSfc  <#- ♦  t^>  i 

In  the  type  of  drying  experiments  reported  here,  L,  the  overt  flow 
path  of  the  vapor  is  a  variable;  it  becomes  longer  as  the  specimen 
"dries-back"  from  the  drying  face. 

If  V  is  the  volume  of  liquid  that  has  been  removed  at  any  time  t, 
and  if  it  is  assumed  that  the  "dry"  end  of  the  specimen  contains  only 
flowing  vapor,  i.e.  the  walls  are  dry  as  the  liquid  front  recedes  into 
the  specimen,  then 

L-Jj  (G8) 

and  then  if  this  value  for  L   is   inserted  in  equation  (G7 )   it  becomes 

dW  .  AV  ,mff_       kBlm       1  (G9) 

dt       ^1/2^  \l  3^      )  v 
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APPENDIX  H 

List  of  Symbols 

A  =  Area  exposed  to  evaporation 

B  ■  Coefficient  of  permeability,  cm  /dyne  -  sec 

C  ■  Concentration  of  gas  at  pressure  P 

C-  «  Constant  for  drying  rate  -  time  curve 

D  ■  Diffusion  coefficient 

E  -  Modulus  of  elasticity 

EF  -  Expansion  factor 

F  ■  Formation  factor 

F  ■  Formation  factor  in  partially  saturated  medium 

G  ■  Geometric  constant  for  diffusion 

H  ■  Intercept  of  drying  time  -  rate  curve 
J(x)  =  Bessel  function 

K  ■  Coefficient  of  permeability,  cm 

K  ■  Absorptivity  coefficient 

8 

K^  ■  Shape  factor  for  diffusion  flow,  1.00 

K  ■  Shape  factor  for  viscous  flow,  2.70 

Le  2 
Kfc  -  Tortuosity,  (■—) 

L  »  Length  of  specimen  in  flow  direction 

L    *  Maximum  permissible  size,  in. 
max 

L  «  Tortuous  length  of  flow  path  through  specimen 
M  «  Molecular  weight 
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N  ■  Slope  of  vacuum  drying  curve 

P  ■  Pressure 

P  ■  Saturation  vapor  pressure 

P  »  Tensile  strength 

P  ■  Average  pressure 

AP  ■  Pressure  difference 

Q  ■  Flow  volume 

R  =  Gas  constant 

R  ■  Resistance  of  the  porous  material  when  saturated  with 
conducting  liquid 

R  ■  Electrical  resistance  of  the  saturant 

S.  ■  Parameter  corresponding  to  fraction  of  impinging  molecules 
undergoing  diffuse  collision  with  wall,  0.75 

S  ■  Degree  of  saturation 

T  ■  Absolute  temperature 

V.  ■  Accommodating  volume 

V  »  Volume  increase  due  to  freezing 

V  ■  Volume  of  pores 
p  t~ 

V_  ■  Volume  of  space  available  to  expelled  water  in  paste 

W  =  Average  moisture  distribution 
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a  sc  Maximum  radius  of  pore 

b  =  Slope  of  drying  time  -  rate  curve 

d  =  Pore  diameter 

QL  m  Rate  of  movement  of  freezing  front 
dt 

e  =  Percent  expansion  at  i  cycles 

h  =  Specified  number  of  cycles,  300 »  at  which  the  freeze-thaw 
exposure  is  to  be  terminated 

i  *  Number  of  cycles  at  which  e  reaches  one-tenth  of  one 
percent  or  300  cycles,  whichever  is  fewer 

m  *=  Hydraulic  radius 

r  =  Pore  radius 

r'  =  Critical  nucleation  radius 

t  =  time 

toe  =  Time  of  drying  at  95  percent  relative  humidity 

tg  =  Time  of  saturation 

u  «=  Apparent  linear  velocity  at  P 

v  =  Mean  thermal  molecular  velocity 

y  =  Ordinate  of  the  pore  size  distribution  curve 

a  =  Positive  root  of  the  Bessel  function  of  the  first  kind 
of  order  zero 

Y  =  Surface  free  energy  (surface  tension) 

«r  =  Porosity,  by  volume 

\  =  Capillary  porosity  of  the  paste 

1   =  Viscosity 

©  =  Constant  angle 

X  =  Cylinder  length 

H-    =  Poisson*s  ratio 

v   =  Function  of  r 
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P    =  Density 

r  =  Thickness  of  the  paste  shell  around  an  aggregate 
particle 

X  =  Rate  of  evaporation 
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APPENDIX  I 

Assumptions  of  Analysis  of  Air-Drying  by  Diffusion 
It  can  be  shown  that  the  series  of  equation  (13)  is  very  closely 
approximated  by  only  the  first  term  when  the  time,  t,  is  sufficiently 
large.  Analysis  of  the  conditions  of  the  air-drying  experiments  during  the 
second  falling  rate  period  shows  that  this  condition  is  more  than  satisfied. 
When  the  first  term  of  equation  (13)  is  differentiated  the  result  is 


£--«« 


[^8  exp  ("  D"2t)J         (ll> 


which  is  equal  to 


|*"»  -  D*2  W  (12) 


The  boundary  conditions  imposed  on  the  solution  of  equation  (13)  and, 
consequently,  on  equation  (12)  are  not  entirely  harmonious  with  the  flow 
process  that  might  be  assumed  to  be  occurring  during  air-drying.  The 
concentration  function,  f(r),  in  equation  (10)  must  be  continuous  between 
r  ■  0  and  r  =■  a.  The  general  picture  of  the  physical  process  {hk,   67)  is 
one  of  a  moving  boundary  between  the  vapor  and  liquid  phases  as  the  liquid 
front  retreats  back  into  the  sample.  Thus  the  term  for  the  length  of 
vapor  flow  path,  a,  is  variable  and  not  a  constant,  as  is  assumed  in 
equation  (13). 

If  this  is  the  process  actually  occurring  during  the  later  stages  of 
air  drying,  the  application  of  equations  (13)  and  (12)  is  inappropriate. 
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The  remainder  of  this  Appendix  is  a  discussion  of  the  reasons  that  seem 
to  justify  the  use  of  these  equations  and,  therefore,  to  negate  the 
simple  moving  boundary  process  just  described. 

Equation  (12)  shows  a  linear  relationship  between  the  rate  of  drying, 
dW/dt,  and  the  moisture  content,  W.   If  the  boundary  were  moving,  one 
would  expect  the  data  to  generate  incrementally  short  segments  of  differing 
slope  so  that  a  non-linear  curve  would  result. 

The  data  from  this  study  do  seem  to  show  a  linear  relationship 
between  the  rate  of  drying  and  the  water  content,  as  indicated  by  equation 
(12).  Fig.  II  shows  an  example  of  typical  data. 

Three  hypotheses  can  be  postulated  as  possible  explanations.  One 
is  that  the  data  are  really  not  linear,  but  only  seem  to  be  on  the  scale 
plotted.   Such  a  possibility  can  only  be  established  with  farther  data 
and  statistical  analysis.   For  the  purposes  here  it  must  be  assumed  that 
the  data  are  what  they  seem,  i.e.  linear.  The  fact  that  the  data  extend 
over  large  differences  in  water  content  support  the  conclusion  of 
linearity. 

Other  workers  have  obtained  data  that  they  consider  to  be  linear  in 
the  variables  in  question.  Hancox  (Vf)  showed  a  linear  relationship 
between  moisture  content  and  air-drying  rates  for  hardened  Portland  cement 
pastes  in  the  second  falling  rate  period  and  analyzed  his  data  by  means  of 
the  one -dimensional  analog  of  equation  (12)  (equation  5)>  which  contains 
the  same  assumptions  discussed  above. 

A  second  possibility  is  a  variable  coefficient  of  diffusion  that  is 
inversely  proportional  to  the  movement  of  the  boundary.  Aside  from  no 
obvious  physical  reason  for  such  behavior,  it  would  be  necessary  to  assume 
that  the  diffusion  coefficient  changes  just  enough  to  keep  the  data  linear, 
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which  would  be  coincidental,  to  say  the  least. 

The  third  possibility,  which  seems  the  most  likely,  is  that  the  flow 
path  length  is  relatively  constant  at  this  rather  advanced  stage  of  the 
drying  process,  or  at  least  that  the  change  in  flow  path  is  small  compared 
with  the  path  length.   One  can  envisage  a  physical  state  in  which  the 
liquid  has  retreated  until  it  occupies  small  pores  and  surface  irregulari- 
ties in  the  void  system.   Subsequent  drying  would  involve  evaporation 
from  these  surfaces  and  a  vapor  flow  by  diffusion  along  a  path  of  sensibly 
constant  length.  The  result  would  be  a  relationship  that  would  satisfy 
the  conditions  leading  to  equation  (12),  and  to  data  such  as  are  in  fact 
found  in  this  study,  and  others  (^7). 

From  the  above  discussion  and  the  findings  of  this  study,  further 
work  to  establish  the  actual  process  occurring  during  the  later  stages 
of  air-drying  is  indicated. 
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